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ABSTRACT 

About one-third of the ~1200 transiting planet candidates detected in the first four months of 
Kepler data are members of multiple candidate systems. There are 115 target stars with two candidate 
transiting planets, 45 with three, 8 with four, and one each with five and six. We characterize the 
dynamical properties of these candidate multi-planet systems. The distribution of observed period 
ratios shows that the vast majority of candidate pairs are neither in nor near low-order mean motion 
resonances. Nonetheless, there are small but statistically significant excesses of candidate pairs both 
in resonance and spaced slightly too far apart to be in resonance, particularly near the 2:1 resonance. 
We find that virtually all candidate systems are stable, as tested by numerical integrations that 
assume a nominal mass-radius relationship. Several considerations strongly suggest that the vast 
majority of these multi-candidate systems are true planetary systems. Using the observed multiplicity 
frequencies, we find that a single population of planetary systems that matches the higher multiplicities 
underpredicts the number of singly-transiting systems. We provide constraints on the true multiplicity 
and mutual inclination distribution of the multi-candidate systems, revealing a population of systems 
with multiple super-Earth-size and Neptune-size planets with low to moderate mutual inclinations. 
Subject headings: Celestial mechanics - planetary systems; Planets and satellites: dynamical evolution 

and stability; Planetary systems; Planets and satellites: fundamental parameters; 

Planets and satellites: general; Facilities: Kepler Space Telescope 



1. INTRODUCTION 

Kepler is a 0.95 m aperture space telescope that 
uses transit photometry to determine the frequency 
and characteris tics o f planets and pla n etary systems 
jBorucki et al. 120101: iKoch et al.l 120101: iJenkins et all 
[2010.: .Caldweh et al.l I2010D . The focus of this NASA 
Discovery mission's design was to search for small plan- 
ets in the habitable zone, but Kepler^s ultra-precise. 
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long-duration photometry is also ideal for detecting sys- 
tems w ith multiple transiting planets. iBorucki et al.l 
(|2011aD presented 5 Kepler targets with multiple tran- 
siting planet candidates, and these candidat e planetary 
systems were analyzed bvl Steffen et a l.' f2010'). A system 
of three plaiiets (K epler-9=KOL377,[Holman et al.' 201(| 
[Torres et ai][201ll) was then reported, in which a near- 
resonant effect on the transit times allowed Kepler data 
to confirm two of the planets and characterize the dy- 
namics of the syster n. A compact system of six transiting 
planets (Kepler-11. iLissauer et aLll2011| ). five of which 
were confirmed by their mutual gravitational interactions 
exhibited through transit timing variations (TTVs), has 
also been reported. These systems provide important 
data for understanding the dynamics, formation and evo- 
lution of_glBneta2;_sj;2teins ([Koch et al. 1996; Agol et aD 
12005': 'Holman & Murrav 2005: Fabrvckvl l2009i : iHolma^ 
EoiO; Ragozzine & Holman 2010). 

A catalog of 1235 candidate planets evident in the 
firs t four and a half mon ths of Kepler data is presented 
in (jBorucki et al.l[2"011bl henceforth Bll). Our analy- 
sis is based on the data presented by Bll, apart from 
the following modifications: Five candidates iden tified 
as very likely false positives bv [Howard "etall([20ll . Ke- 
pler Objects of Interest (KOIs) 1187.01, 1227.01, 1387.01, 
1391.01 and 1465.01, are removed from the list. The ra- 
dius of KOI-1426.03 is reduced from 35 R® to 13 R® 
(Section 4). Bll reported that KOI-730.03's period was 
almost identical to that of KOI-730.02, putting these 
two planets in the 1:1 resonance. Upon more careful 
lightcurve fits, we disfavor the co-orbital period of KOI- 
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730.03 near 10 days as an alias. The period of KOI-730.03 
is therefore increased to 19.72175 and its radius increased 
to 2.4 R^. We do not consider the 18 candidates for 
which only a single transit was observed (and thus the 
period is poorly-constrained), nor the 13 remaining can- 
didates with estimated radii Rp > 22 A R^. None of 
the removed candidates is in a multi-candidate system, 
but both of the candidates with adjusted parameters are. 
We are left with a total of 1199 planetary candidates, 408 
of which are in multiple candidate systems, around 961 
target stars. 

Bll also summarizes the observational data and follow- 
up programs being conducted to confirm and character- 
ize planetary candidates. Each object in this catalog is 
assigned a vetting flag (reliability score), with 1 signify- 
ing a confirmed or validated planet (>98% confidence), 
2 being a strong candidate (>80% confidence), 3 being 
a somewhat weaker candidate (>60% confidence), and 
4 signifying a candidate that has not received ground- 
based follow-up (this sample also has an expected relia- 
bility of >60%). For most of the remainder of the paper, 
we refer to all of these objects as "planets", although 
99% remain unvalidated and thus should more properly 
be terr ned "planet ca ndidates." iRagozzine & Holmad 
(|2010( ). iLissauer et"all (|2011l ). andEatham ct al. (2011^ 
point out that KOIs of target stars with multiple planet 
candidates are less likely to be caused by astrophysical 
false positives than are KOIs of stars with single candi- 
dates, although the possibility that one of the candidate 
signals is due to, e.g., a background eclipsing binary, can 
be non-negligible. 

The apparently single systems reported by Bll are 
likely to include many systems that possess additional 
planets that have thus far escaped detection by being too 
small, too long-peri od, or too h i ghly in clined with respect 
to the line of sight. iFord et all (|20lH ) discuss candidates 
exhibiting transit timing variations (TTVs) that may be 
caused by these non-transiting planets, but so far no non- 
transi ting planets have been identified. iLatham et al.l 
compare various orbital and physical properties 
of planets in single candidate systems to those of plan- 
ets in multi-candidate systems. We examine herein the 
dynamical aspects of Kepler^s multi-planet systems. 

Multi-transiting systems provide numerous insights 
that are difficult or impossible to gain from single- 
transiting systems (jRagozzine &: HolmanI [2010l ) . These 
systems harness the power of radius measurements from 
transit photometry in combination with the illuminat- 
ing properties of multi-planet orbital architecture. Ob- 
servable interactions between planets in these systems, 
seen in TTVs, are much easier to characterize, and our 
knowledge of both stellar and planetary parameters is 
improved. The true mutual inclinations between planets 
in multi-transiting systems, while not directly observ- 
able from the light curves themselves, are much easier 
to obtain than in other syste ms, both individuall y and 
statistically (see SectionlSland lLissauer et "aLll2011[ ). The 
distributions of physical and orbital parameters of plan- 
ets in these s ystems are invaluab le for comparative plan- 
etology (e.g.. iHavel et al.ll2011| ) and provide important 
insights into the processes of planet formation and evo- 
lution. It is thus very exciting that Bll report 115 dou- 
bly, 45 triply, 8 quadruply, 1 quintuply and 1 sextuply 
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Characteristics of Systems with Two Transiting Planets 
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Fig. 1. — Planet period vs. radius for all 1199 planetary candidates from Bll that we consider herein. Those planets that are the only 
candidate for their given star are represented by black dots, those in two-planet systems as blue circles (open for the inner planets, filled 
for the outer ones), those in three-planet systems as red triangles (open for the inner planets, filled for the middle ones, filled with black 
borders for the outer ones), those in four planet systems as purple squares (inner and outer members filled with black borders, second 
members open, third filled), the five candidates of KOI-500 as orange pentagons and the six planets orbiting KOI-157 (Kepler-11) as green 
hexagons. It is immediately apparent that there is a paucity of giant planets in multi-planet systems; this difference in the size distributions 
is quantified and discussed by Latham ct al. (2011). The upward slope in the lower envelope of these points is caused by the low SNR of 
small transiting planets with long orbital periods (for which few transits have thus far been observed). Figure provided by Samuel Quinn. 



We begin by summarizing the characteristics of ife- 
pZer's muhi-planet systems in Section 2. The reliabihty 
of the data set is discussed in Section 3. Sections 4-6 
present our primary statistical results. Section 4 presents 
an analysis of the long-term dynamical stability of can- 
didate planetary systems for nominal estimates of plan- 
etary masses and unobserved orbital parameters. Evi- 
dence of orbital commensurabilities (mean motion reso- 
nances) among planet candidates is analyzed in Section 
5; this section also includes specific comments on some of 
the most interesting cases in the reported systems. The 
characteristic distributions of numbers of planets per sys- 
tem and mutual inclinations of planetary orbits are dis- 
cussed i n Section 6. We compare the Kepler-11 planetary 
system (jLissauer et al.ll2011 ) to other Kepler candidates 
in Section 7. We conclude by discussing the implications 
of our results and summarize the most important points 
presented herein. 

2. CHARACTERISTICS OF KEPLER MULTI-PLANET 
SYSTEMS 

The light curves derived from Kepler data can be used 
to measure the radii and orbital periods of planets. Pe- 
riods are measured to high accuracy. The ratios of plan- 
etary radii to stellar radii are also well measured, ex- 
cept for low signal to noise transits (e.g., planets much 
smaller than their stars, those for which few transits are 
observed, and planets orbiting faint or variable stars). 
The primary source of errors in planetary radii are prob- 
ably errors in estimated stellar radii, which are calcu- 
lated using photometry from the Kepler Input Catalog 
(Bll). In some cases, the planet's radius is underesti- 
mated because a significant fraction of the flux in the 
target aperture is provided by a star other than the one 
being transited by the planet. This extra light can of- 



ten be dete rmined and correct ed for using data from the 
spacecraft (jBrvson et al.|[2010l ) and/or from the ground, 
although dilut ion by faint physic ally-bound stars is diffi- 
cult to detect (jTorres et al.ll20li[ ). Spectroscopic studies 
can, and in the future will, be used to estimate the stel- 
lar radii, thereby allowing for improved estimates of the 
planetary radii; this procedure will also provide more ac- 
curate estimates of the stellar masses. 

We list the key observed properties of Kepler two- 
planet systems in Table [1] three-planet systems in Ta- 
ble El four-planet systems in Table |3l and those systems 
with more than four planets in Table HI Planet indices 
given in these tables signify orbit order, with 1 being the 
planet with the shortest period. These do not always cor- 
respond to the post-decimal point portion (.01, .02, etc.) 
of the KOI number designation of these candidates given 
in Bll, for which the numbers signify the order in which 
the candidates were identified. These tables are laid out 
differently from one another because of the difference in 
parameters that are important for systems with differ- 
ing numbers of planets. In addition to directly observed 
properties, these tables contain results of the dynamical 
analyses discussed below. 

We plot the period and radius of each active Ke- 
pler planetary candidate observed during the first four 
and one-half months of spacecraft operations, with spe- 
cial emphasis on multi-planet systems, in Figure [TJ We 
present galleries of multi-planet candidates, representing 
the planetary sizes and periods, in Figures [JHH The 
ratio(s) of orbital periods is a key factor in planetary 
system dynamics; the cumulative distributions of these 
period ratios for various classes of planetary pairings are 
plotted in Figure [5j Figure |6] compares the cumulative 
distribution of the period ratio of neighboring pairs of 
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TABLE 2 

Characteristics of Systems with Three Transiting Planets 
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2.785831 


26.8 
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1.59 


1.77 
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54.9 
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2.108971 


27.9 
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3.17 
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23.9 
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18.5 


3.29 


2.563540 


17.3 


806 


3.10 


12.05 


2.068502 


8.3 


8.99 


2.373352 


8.7 


812 


2.45 


2.19 


6.005814 


39.0 


2.12 


2.302249 


20.8 


829 


1.93 


2.63 


1.912284 


19.4 


2.72 


2.067646 


19.5 


864 


2.17 


1.82 


2.265275 


25.5 


1.79 


2.052797 


24.3 


884 


1.24 


3.00 


2.829393 


27.8 


2.74 


2.169267 


18.1 


898 


2.41 


3.03 


1.889897 


13.8 


2.54 


2.056092 


15.3 


899 


1.28 


1.69 


2.151439 


22.7 


1.72 


2.160301 


20.8 


907 


2.07 


3.52 


3.446956 


30.7 


3.41 


1.824674 


13.7 


921 


1.48 


2.30 


2.717125 


31.1 


2.68 


1.762252 


15.2 


934 


3.16 


2.02 


2.130198 


20.2 


2.41 


1.510392 


12.6 


935 


3.59 


3.16 


2.043775 


16.5 


2.47 


2.055839 


18.8 


941 


3.43 


5.39 


2.762270 


18.3 


6.61 


3.747567 


19.1 


961 


14.44 


3.94 


2.677706 


8.5 


10.72 


1.536625 


4.6 


1306 


2.14 


2.20 


1.930557 


20.5 


1.81 


1.705578 


17.6 


1422 


2.02 


3.06 


1.613089 


10.5 


2.92 


3.398094 


23.4 


1426 


3.29 


7.74 


1.926989 


10.3 


12.30 


2.002740 


7.5 



TABLE 3 

Characteristics of Systems with Four Transiting Planets 



KOI # 


Rp.i (R®) 


«P,2 (R®) 


P2/P1 


Ai,2 


Rp.3 (R®) 


P3/P2 


A2,3 


RvA (R®) 


Pi/Ps 


A3, 4 


70 


1.60 


0.60 


1.649977 


22.6 


2.27 


1.779783 


20.8 


1.97 


7.150242 


53.5 


117 


1.29 


1.32 


1.541373 


19.4 


0.68 


1.623476 


25.1 


2.39 


1.853594 


22.3 


191 


1.39 


2.81 


3.412824 


35.8 


11.56 


6.350795 


20.2 


1.49 


1.258207 


2.9 


707 


2.19 


3.36 


1.652824 


13.4 


2.48 


1.459577 


9.8 


2.64 


1.290927 


7.3 


730 


1.83 


2.29 


1.333812 


9.4 


3.09 


1.500997 


11.0 


2.63 


1.333948 


7.5 


834 


1.39 


1.39 


2.943873 


44.2 


1.85 


2.149832 


28.5 


4.87 


1.787503 


12.5 


880 


2.01 


2.76 


2.476866 


25.6 


4.90 


4.480164 


28.7 


5.76 


1.948715 


11.1 


952 


1.14 


2.27 


2.037706 


19.5 


2.32 


1.483153 


9.2 


2.44 


2.602734 


21.3 



Kepler planet candidates with the comparable distribu- 
tion for radial velocity planets. The derivatives of these 
distributions with respect to period ratio, which show 
spikes near common period ratios, are displayed in Fig- 
ure [71 

The probability that a Kepler target star hosts at 
least one detected transiting planet is 961/160171 ~ 
0.006, while the probability that a star with one de- 
tected transiting planet hosts at least one more is much 
higher, 170/961 ~ 0.177. To determine whether or not 
the observed transiting planets are randomly distributed 
among host stars, we compare the observed distribution 
of planets per system to three random populations. In 



the first of these synthetic populations, each of 1199 plan- 
ets is randomly assigned to one of 160171 stars. The 
number of stars containing j planets in this random pop- 
ulation is equivalent to a Poisson distribution with a 
mean of 1199/160171 « 0.0075, and is denoted as "Pois- 
son 0" . For the second synthetic population, "Poisson 1" , 
we again use a Poisson distribution, also forced to match 
the observed total number of Kepler planets, 1199, but 
in this case the second requirement is to match the to- 
tal number of planetary systems (stars with at least 1 
observed planet), 961, rather than the total number of 
stars, 160171. Results are shown in Tabled which also 
includes comparisons to a third random distribution dis- 
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TABLE 4 

Characteristics of Systems with Five or Six Transiting Planets 



property 


KOI-500 


KOI-157 (nominal Koplcr-11) 


Kcplcr-11 (Lissauer et al. 2011) 


Rp.i (Rffi) 


1.23 


1.70 


1.97±0.19 


Rp,2 (R®) 


1.47 


3.03 


3.15±0.30 


P2/P1 


3.113327 


1.264079 


1.26410 


A12 


40.5 


6.8 


7.0 


Rp,3 (Re) 


2.06 


3.50 


3.43±0.32 


P3/P2 


1.512077 


1.741808 


1.74182 


A23 


12.7 


13.1 


15.9 


(Re) 


2.75 


4.21 


4.52±0.43 


P4/P3 


1.518394 


1.410300 


1.41031 


A34 


10.4 


7.3 


10.9 


Rp.5 (Re) 


2.83 


2.22 


2.61±0.25 


P5/-P4 


1.349929 


1.459232 


1.45921 


A45 


6.8 


8.8 


13.3 


flp.e (Re) 




3.24 


3.66±0.35 






2.535456 


2.53547 


A56 




24.1 
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Fig. 2. — Gallery of candidate planetary systems of four to six 
planets. Each horizontal line represents a separate system, as la- 
beled. They are sorted by, first, the number of candidates, and 
second, the innermost planet's orbital period. The size of each 
dot is proportional to the size of the planet that it represents, and 
in each system the size orderings move from hot (large) to cool 
(small) colors (red for the largest planet within its system, then 
gold, green, aqua, and, if additional planets are present, navy and 
gray). There is a clear trend for smaller planets to be interior to 
larger planets, but this is due to the greater detectability of small 
planets at shorter orbital period. 

cussed below and to the population of planets detected by 
radial velocity surveys normalized to the same total num- 
ber of planets as the Kepler candidates considered here. 
There are sharp differences between Kepler observations 
and Poisson 0: The number of observed Kepler multi- 
transiting systems greatly exceeds the random distribu- 
tion, showing that transiting planets tend to come in 
planetary systems, i.e., once a planet is detected to tran- 
sit a given star, that star is much more likely to be or- 
bited by an additional transiting planet than would be 
the case if transiting planets were randomly distributed 
among targets. Poisson 1, the random distribution that 
fits the total number of planets and the total number 
of stars with planets, underpredicts the numbers of sys- 
tems with three or more planets, suggesting that there is 
more than one population of planetary systems. To carry 
the study further, wc used a third random distribution, 
"Poisson 2", constrained to fit the observed numbers of 
multi-planet systems (170) and planets within said sys- 
tems (408). The Poisson 2 random distribution fits the 



numbers of Kepler systems with multiple transiting plan- 
ets quite well, but it only accounts for ~ 35% of the single 
planet systems detected. As Kepler does not detect all 
planets orbiting target stars for which some planets are 
observed, the observed multiplicities differ from the true 
multiplicity; we explore these differences in Section [51 

ft is difhcult to estimate the corresponding distribution 
for planets detected through Doppler radial velocity (RV) 
surveys, since we do not know how many stars were spec- 
troscopically surveyed for planets for which none were 
found. Ho wever, a search th rough the Exoplanet Orbit 
Database (|Wright et al.|[20Tll) as of January 2011 shows 
that 17% of the planetary systems detected by RV sur- 
veys are multi-planet systems (have at least two plan- 
ets), compared to 21% for Kepler. Including a linear RV 
trend as evide nce for an additional planetary companion, 
IWright et all (|20093 show that at least 28% of known RV 
planetary systems are multiple. Note, however, that the 
radial velocity surveys are generally probing a population 
of more massive planets and one that includes longer 
period planets than the Kepler sample discussed here. 
Some of the most recent RV planets are from the same 
population of small-pl anet, short-period mul t iples seen in 
the K epler data (e.g., iLo Curto et al.|[2010t iLovis et all 
120111 ). Also, the observed multiplicity of RV systems is 
far less sensitive to the mutual inclinations of planetary 
orbits than is the multiplicity of Kepler candidates. The 
final column in Table [S] shows the distribution of radial 
velocity detected planets normalized to the total number 
of Kepler candidates. 

In systems with multiple transiting candidates, the 
ratio(s) of planetary periods and the ratio(s) of plane- 
tary radii arc well-measured independently of uncertain- 
ties in stellar properties. This suggests an investigation 
into whether outer planets are larger than inner plan- 
ets on average, in order to provide constraints on the- 
ories of planet formation. (The ratio of orbital periods 
is discussed in Section [5l) When considering the entire 
multi-candidate population, there is a slight but signifi- 
cant preference for outer planets to be larger (Figure [T]). 
However, since planets with longer periods transit less 
frequently, all else being equal the SNR scales as P^^/"^, 
after accounting for the increase in SNR due to longer du- 
ration transits (assuming circular orbits). To debias the 
radius ratio distribution, each planetary system is inves- 
tigated and if the smallest planet cannot be detected at 



6 



Lissauer et al. 



TABLE 5 



# planets per system 


Kepler 


Poisson 


Poisson 1 


Poisson 2 


RV (sealed to Kepler) 





159210 


158976.5 


1630.9 


227.3 




1 


791 


1190.1 


755.5 


231.1 


836.8 


2 


115 


4.5 


175.0 


117.6 


97.2 


3 


45 


1.1x10"^ 


27.0 


39.9 


27.8 


4 


8 


2.1x10"'°'' 


3.1 


10.1 


10.4 


5 


1 


3.1x10"** 


2.9x10"^ 


2.1 


3.5 


6 or more 


1 


3.9x10"" 


2.4x10"^ 


4.1x10"^ 


3.5 


Total stars: 


160171 


160171 


2591.9 


628.4 




Total planets: 


1199 


1199 


1199 


639.1 


1199 


Total stars with planets; 


961 


1194.5 


961 
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981.0 
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Fig. 3. — Three planet candidate systems; same format as Fig- 
ure [2] 



SNR > 16 if it were placed in the longest period, that 
planet is removed. The remaining population of systems 
are not sensitive to observational bias, since each planet 
could be detected at all periods. The distribution of ra- 
dius ratios (i?p^o/^p,i), where the subscripts o and i refer 
to the inner and outer members of the pair of planets, for 
debiased planetary systems is shown by the solid curve 
in Figure |8l Interestingly, neighboring planets tend to 
have very similar radii, with most of the population near 
Rp^o/Rp^i ~ 1. This strong tendency is illustrated by 
comparing the cumulative distribution of observed radii 
ratios with ratios of radii randomly drawn from the debi- 
ased distribution (dashed curve in Figure H]). Since there 
is no requirement that the apparent depths of false posi- 
tives and/or planets to be comparable, this suggests that 
false positives are not common among systems with radii 
ratios near unity. In the debiased distribution, there is 
no significant preference for the outer planet to be larger 
than the inner one. 

Examination of this debiased sample for additional 
non-trivial trends and correlations (between radius ratio, 
period ratio, radius, period, multiplicity, near-resonance 
(defined by (i discussed in Section 5), and stellar proper- 
ties, reveals a few additional clues into the distribution of 
the Kepler candidates. When there is a giant {R > 6Rq) 
planet, it is usually, but not always, at longer period. Ra- 
dius ratios near 1 are almost exclusively found for plan- 
ets less than 4 i?^. Another trend is that whenever the 
period ratio between two neighboring detectable plan- 
ets exceeds 5, the sizes between the planets were similar 
and, conversely, large radius ratios only occurred for pe- 
riod ratios <3. Furthermore, large period ratio systems 
are only found around small stars (i?^ < Rq). These 
latter examples may be due to observational bias. 

3. RELIABILITY OF THE SAMPLE 

Very few of the planetary candidates presented herein 
have been validated or co nfirme d as true exoplanets. As 
discus sed in Bll. Brownl (|2003[ ) and iMorton fc JohnsonI 
(|2011f ). the overwhelming majority of false positives are 
expected to come from eclipsing binary stars. Many of 
these will be eclipses of fainter stars (physically associ- 
ated to the target or background/faint foreground ob- 
jects) within the aperture, others grazing eclipses of the 
target star. Similarly, transits of fainter stars that ap- 
pear within the target aperture by large planets may also 
mimic small planets transiting the target star. The vast 
majority of all active Kepler planet candidates were iden- 
tified using automated programs that do not discriminate 
between single candidates and systems of multiple plan- 
ets. In addition, < 1% of single and multiple candidate 
systems were identified in subjective visual searches; this 
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Fig. 4. — Two planet candidate systems; same format as Figure [J] 
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Fig. 5. — Cumulative number, N, of pairs of Kepler planets or- 
biting the same star with period ratio, P (= Po/Pi), less than 
the value specified. Black dashed curve shows all pairings, solid 
black curve shows pairings of neighboring planets; solid red curve 
shows all pairings in two candidate systems; dashed blue curve 
shows all pairings in three-planet systems, solid blue curve shows 
all adjacent pairings in three-planet systems, solid light green curve 
shows inner pair of planets in three-planet systems, the solid dark 
green curve shows outer pair of planets in three-planet systems, the 
dashed pink curve shows all pairings in the 4, 5 and 6 planet sys- 
tems, solid pink curve shows neighboring pairings in the 4, 5 and 
6 planet system. Sixty-five pairs of planets in the same system, 
including 31 adjacent pairs, have P > 5 and are not represented 
on this plot. 



Fig. 7. — Slope of the cumulative fraction of Kepler neighboring 
planet pairs (solid black curve) and multi-planet systems detected 
via radial velocity (dashed red curve) with period ratio less than 
the value specified. The slope for the Kepler curve was computed 
by taking the difference in period ratio between points with 
differing by 4 and dividing this difference by 4. The slope for the 
RV curve was computed hy taking the difference in period ratio 
between points with N differing by 3 and dividing by 3, and then 
normalizing by multiplying the value by the ratio of the number of 
Kepler pairings to the number of radial velocity pairings (3.9). The 
spikes in both curves near V = 2 and the Kepler curve near V = 
1.5 show excess planets piling up near period commensurabilities. 
Two points on the Kepler curve near period ratio 1.5 lie above 
the plot, with the highest value being 299.84; this sharp peak is 
produced by the excess of observed planet pairs in or very near the 
3:2 mean motion resonance. 




1 10 100 



Period Ratio 

Fig. 6. — Cumulative fraction of neighboring planet pairs for 
Kepler candidate multi-planet systems with period ratio less than 
the value specified (solid curve) . The cumulative fraction for neigh- 
boring pairs in multi-planet systems detected via radial velocity is 
also shown (dashed curve), and includes data from exoplanets.org 
as of 29 January 2011. (a) Linear horizontal axis, as in Figure [5] 
The data are normalized for the number of adjacent pairs with V 
< 5, which equals 207 for the Kepler candidates and 28 for RV 
planets, (b) Logarithmic horizontal axis. All 238 Kepler pairs are 
shown; the three RV pairs with V > 200 are omitted from the plot, 
but used for the calculation of Ntot, which is equal to 61. 




Fig. 8. — Cumulative distribution of the ratio of planetary radii 
{Rp,o/ Rp,i) for neighboring pairs of transiting planet candidates 
for which both planets are detectable at the longer period. The 
distribution of radii ratios for all such neighboring pairs (solid line, 
71 ratios) is the same (within statistical uncertainty) as that of 
such pairs of planet candidates that are near the 2:1 mean motion 
resonance (dotted line, 22 ratios). Interestingly, neighboring plan- 
ets tend to have very similar radii, with most of the population 
near Rp.o/Rp,i ~ 1. This strong tendency is illustrated by show- 
ing for contrast the cumulative distribution of ratios between radii 
randomly drawn from the debiased distribution (dashed line). In 
the debiased distribution, there is no significant preference for -Rp,o 
to be greater or less than Rp^i- 
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fraction is too small to significantly affect our statistical 
results. 

The fact that planetary candidates are clustered 
around targets — there are far more targets with more 
than one candidate than would be expected for randomly 
distributed candidates (Section [5]) — suggests that the 
reliability of the multi-candidate systems is likely to be 
higher than that of the single candidates. We know from 
radial velocity obser vations that planets frequently come 
in multiple systems ([Wright et al.l 120091 ) . whereas astro- 
physical false positives are expected to be nearly random. 
The mere presence of multiple candidates increases our 
confidence that most or all are real planets, since the 
probability of multiple false positive signals is the prod- 
uct of the probabilities of two o r more relatively rare 
cases (IRagozzine fc Holmanll2010l ). (Triple star systems 
would be dynamically unstable for most of the orbital pe- 
riod ratios observed. Their dynamics would also give rise 
to ve ry large eclipse timing variations fe.g..lCart er et ahl 
1201 It ), which are generally not seen.) Additionally, the 
concentration of candidate pairings with period ratio 
near first-order mean motion resonances such as 3:2 and 
2:1 (Figure [S] and [T]) suggests that these subsamples likely 
have an even larger fraction of true planets, since such 
concentrations would not be seen for random eclipsing 
binaries. These qualitative factors have not yet been 
quantified. 

Strong TTV signals are no more common for planets in 
observe d multi-transiting candidates than for single can- 
didates (jFord et al.ll201lD . and it is known that some of 
these systems are actually stellar triples masquerading as 
planetary systems, suggesting that large TTVs are often 
the result of stellar systems ma squerading as planetary 
systems fe.g.. lSteffen et al.1l2011[ ). Nonetheless, it is pos- 
sible that weaker TTVs correlate with multi-transiting 
systems, but that these sm all TTVs (e.g.. thos e in the 
Kepler-11 planetary svstem. fLissauer et al.|[20Tl[) require 
more data and very careful study to reveal. 

The ratio of transit durations normalized to the ob- 
served periods of candidates in multi-transiting systems 
can also be used to identify false positives (generated by a 
blend of two objects eclipsing two stars of different den- 
sities) or as a probe of the eccentricity and inclination 
distributions. Transit durations depend upon orbital pe- 
riod, and it is convenient to normalize values for candi- 
dates in a system by dividing measured duration by P^^^; 
this normalization would yield the same value for plan- 
ets of differing periods on circular orbits with the same 
impact parameter around the same star. For this study, 
we assumed circular orbi ts and a random dist ribution of 
impact parameters as in iHolman et al.l ()2010( ) . We find 
that the distribution of normalized transit duration ra- 
tios is consistent with that of a population of planetary 
systems with no contamination. The consistency of nor- 
malized duration ratios is not a strong constraint for false 
positives, but the observed systems also pass this test. 

There are five pairs for which the observed transit du- 
ration ratio falls in the upper or lower 5% of the distri- 
bution predicted. These are: KOI-864.01 and 864.03, 
291.02 and 291.01, 1426.01 and 1426.03, 645.01 and 
645.02 and 1089.02 and 1089.01; independent evidence 
from the shape of the transit light curve indicates that 
KOI-1426.03's transits are grazing (Section|4|). However, 



this test is not definitive enough to rule out the plan- 
etary interpretation, in favor of a blended eclipsing bi- 
nary hypothesis, for any of these pairs. Given 238 pairs 
of neighboring planets among the multiple planet candi- 
date systems identified by Kepler, a 10% false alarm rate 
(corresponding to the extreme 5% on either side of the 
distribution) could easily result in nearly two dozen such 
pairs of planets simply due to chance. The substantially 
smaller fraction of extreme normalized duration ratios 
that are observed may result from the lower SNR (and 
thus decreased likelihood for detection) of grazing and 
near-grazing transits. In sum, we do not find evidence 
for false positives based on the ratio of transit durations 
of KOIs with a common host star. 

4. LONG-TERM STABILITY OF PLANETARY SYSTEMS 

Kepler measures planetary sizes, whereas masses are 
the key parameters for dynamical studies. For our 
dynamical studies, we convert planetary radii, Rp, to 
masses, Mp, using the following simple formula: 



^Rp_\ 



2.06 



(1) 



where R® and M0 are the radius and mass of the Earth, 
respectively. The power-law of Equation ([T]) was ob- 
tained by fitting to Earth and Saturn; it slightly over- 
estimates the mass of Uranus (17.2 vs. 14.5 M^) 
and slightly underestimates the mass of Neptune (16.2 
vs. 17.1 M0). (Note that Uranus is larger, but 
Neptune is more massive, so fitting both with a mono- 
tonic mass-radius relatio n is not possible.) O bservations 
of transiting exoplanets (jLissauer et al.l [201 IL and refer- 
ences therein) show more significant deviations from the 
relationship given by Equation (1) than do Uranus and 
Neptune, with both denser and less dense planets known, 
but on average the known exoplanets smaller than Saturn 
are consistent with this general trend. The stellar mass, 
taken from Bll, was derived from the R^, and log g esti- 
mates from color ph otometry of the Kepler Input Catalog 
(jBrownet al.l[20ll . 

A convenient metric for the dynamical proximity of the 
j*'* and {j -1-1)*'' planets is the separation of their orbital 
semi-major axes, Oj+i — Oi , measured in units of their 
mutu al Hill sphere radius (|Hil]| 118781 ISmith fc Lissaueil 
[20T0h : 



R 



Hi 



Mj + Mj+i 1 1/3 (aj + Oj+i) 



(2) 



where the index j — (1, Np — 1), with Np being the num- 
ber of planets in the system under consideration, Alj and 
Oj are the planetary masses and semi-major axes, and M^, 
is the central star's mass (the masses of interior planets 
are likely to be much smaller than the uncertainty in the 
star's mass, and thus have been neglected in Equation 

The dynamics of two-planet systems are a special case 
of the 3-body problem that is amenable to analytic treat- 
ment and simple numerically-derived scaling formulas. 
For instance, a pair of planets initially on coplanar cir- 
cular orbits with dynamical orbital separation 



R 



> 2%/3 w 3.46 
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can never de velop crossing; o rbits, and are thus called 
"Hill stable" fGladmanI 119931) . Dynamical orbital sepa- 
rations of all planet pairs in two-planet systems for stellar 
masses given in Bll and planetary masses given by Equa- 
tion ([1]) are listed in Table [TJ We see that all two-planet 
systems obey the criterion given in Expression ([3]) . This 
is additional evidence in favor of a low false positive rate, 
as random false positives would occasionally break this 
criterion. 

The actual dynamical stability of planetary systems 
also depends on the eccen tricities and mutual inclination 
(|Veras fc Armitagd l200l . none of which can be mea- 
sured well from the transit data alone (cf. iFord et aD 
120111) . Circular coplanar orbits have the lowest angu- 
lar mo mentum defici t (AMD), and thus are the most 
stable (jLaplacd 117841 : iLaskaii [19971) . An exception ex- 
ists for planet pairs protected by mean motion reso- 
nances (Subsection 15. 4|) that produce libration of the 
relative longitudes of various combinations of orbital el- 
ements and prevent close approache s. A second ex- 
cepti on exists for retrograde orbits (jSmith fc LissaueJ 
I2009D. but these are implausibl e on cosmogonic grounds 
([Lissauer fc Slartibartfastll2008[ ) . Hence, the Kepler data 
can only be used to suggest that a pair of planets is 
not stable. Instability is particularly likely if their semi- 
ma.jor axes are closer than th e resonance overlap limit 
(|Wisdomlll980t IMalhotralHOQSl) : 



Oi + tto 



2/7 



(4) 



Equation (4) was derived only for one massive planet; 
the other is massless. Conversely, if it is assumed that 
both candidates are true planets, limits within the pa- 
rameter space of planetary masses, eccentricities and mu- 
tual inclinations can be inferred from the requirement 
that the planetary systems are presumably stable for 
the age of the system, which is generally of order 10^^ 
times longer than the orbital timescale for the short- 
period planets that represent the bulk of Kepler's candi- 
dates. However, we may turn the observational problem 
around: The ensemble of duration measurem ents can be 
used to address the eccentricity distribution (iFord et al.l 
[2008); since multi-planet systems have eccentricities con- 
strained by stability re quirements, they can b e used as a 
check to those results ([Moorhead et al.ir2Qll[ ). 

Systems with more than two planets have additional 
dynamical complexity, but are very unlikely to be sta- 
ble unless each neighboring pair of planets satisfies Ex- 
pression ([3]). Nominal dynamical separations (A's) be- 
tween neighboring planet pairs in three-planet systems 
are listed in Table [21 the separations in four-planet sys- 
tems are listed in Table [3l and those for the five-planet 
and six-planet systems in Table [H We find that in the 
overwhelming majority of these cases, the inequality in 
E xpression ([3l) is satisfied by quite a wide factor. 

iSmith fc Lissauei) (|2009D conducted suites of numerical 
integrations to examine the stability of systems consist- 
ing of equal mass planets that were equally spaced in 
terms of A. They demonstrated dynamical survival of 
systems with three comparably-spaced 1 M0 planets or- 
biting a 1 Mq star for 10^" orbits of the inner planet 
in cases where the relative spacing between orbital semi- 
major axes exceeded a critical number (Acrit ~ 7) of 
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Fig. 9. — The orbital separations, expressed in mutual Hill radii 
(Equation O , for the inner (Aj) and outer (Ao) pairs within 3- 
planet systems and adjacent 3-planet sub-systems of the 4-, 5-, and 
6-planet systems of our sample. The colored symbols are the same 
as in Figure [1] and we also include the spacings of planets within 
the Solar System as black circles. Gray region: stability boundary 
for 2-planet systems, A > 2\/3 (Expression |3)l . Dashed line; sta- 
bility boundary used for the simulated population of multi-planet 
systems: A; -|- Ao > 18; some of the nominal systems survived 
long-term integrations despite transgressing that boundary. The 
unstable nominal systems KOI- 191 and KOI-248 both lie within 
the gray region, and all the long-term survivors are outside of it. 

mutual Hill spheres. For systems with five comparably- 
spaced 1 M0 planets, they found Acrit ~ 9. Their calcu- 
lations give slightly larger Acrit for systems of five 0.33 
planets, suggesting a slightly smaller Acrit for the 
somewhat larger planets that dominate the Kepler sam- 
ple of interest here. Systems with planets of substan- 
tially differing masses tend to be less stable because a 
smaller fraction of the AMD needs to be placed in the 
lowest mass planets for orbits to cross. Taking all of 
these factors, as well as the likely presence of undetected 
planets in many of the systems, into account, we consider 
Acrit ~ 9 to be a rough estimate of how close planetary 
orbits can be to have a reasonable likelihood of survival 
on Gyr time scales. 

For each adjacent set of 3 planets in our sample, we 
plot the A separations of both the inner and the outer 
pairs in Figure [S] We find that in some cases A < 9, but 
the other separation is quite a bit larger than Acrit- In 
our simulations of a population of systems in Section [HI 
we impose the stability boundary: 



A,+i,, + A 



> 18, 



(5) 



which accounts for this possibility, in addition to re- 
quiring that inequality in Expression ([3|) is satisfied. 
The form of Expression ([ 5 ] ) is motivated by observations 
within our Solar Svstem (|LissaueHll995[) . 

We investigated long-term stability of all 55 systems 
with three or more planets u sing the hybrid integrator 
within the Mercury package ([Chamberslll999n . We set 
the switchover at 3 Hill radii, but in practice we aborted 
simulations that violated this limit, so for the bulk of the 
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simula tion the rt-body mapping of iWisdom fc Holmad 
(|1991f ) was used, with a time step of 0.05 times as large as 
the orbital perio d of the innermo st planet. The simplest 
implementation (|Nobili fc Roxbu rgh 1986) of general rel- 
ativistic precession was used, an additional potential 



where G is Newton's constant, c is the speed of 
light, and r is the instantaneous distance from the 
star. More sophisticated tr eatments of general relativ- 
ity (jSaha fc Tremainj I1994D are not yet required, due 
to the uncertainties of the masses of the planets and 
stars whose dynamics are being modeled. We neglected 
precession due to tides on or rotational flattening in 
the planets, which are only significant compared to gen- 
eral relativity f or Jupiter-size planets in few-day orbits 
(|Ragozzine fc W oU 2009), and precession due to the ro- 
tational oblateness of the star, which can be significant 
for very close-in planets of any mass. Precession due to 
the time- variable flattening of the star c an generate a sec- 
ular r esonance, compromising stability ()Nagasawa &: LinI 
HOOl); we neglected this effect, treating the star as a point 
mass. We also neglect tidal damping of eccentricities, 
which can act to stabilize systems over long time scales, 
and tidal evolution of semi-major axes, which sometimes 
has a destabilizing effect. We assumed initially circu- 
lar and coplanar orbits that matched the observed peri- 
ods and phases, and chose the nominal masses of Equa- 
tion H]). 

For the triples, quadruples and the five-planet system 
KOI-500, we ran these integrations for 10^" orbital peri- 
ods of the innermost planet, and found the nominal sys- 
tem to be stable for this span in nearly all cases; the two 
exceptions are described below. The most populous tran- 
siting planetary system discovered so far is th e six-planet 
system Kepler-11. In the discovery paper, Lis sauer et al.l 
([201 If) reported a circular model, with transit time vari- 
ations yielding an estimate of the masses of the five in- 
ner planets. They also proposed two more models with 
planets on moderate eccentricities, and slightly different 
masses, which also fit the transit timing signals. Con- 
tinued integration showed that both of those eccentric 
solutions are nominally unstable, at 610 Myr for the all- 
eccentric fit, and at 427 Myr for an integration that was 
a very slight offset from the b/c-eccentric fit. The all- 
circular fit remained stable for the duration of our 1 Gyr 
integration, so we conclude that the fits of that paper 
are physically plausible. Introducing a small amount of 
eccentricity can compromise stability, but our results do 
not exclude eccentricities of a few percent because (1) 
planetary masses have significant uncertainties, and we 
have not conducted a survey of allowed mass/eccentricity 
combinations, and (2) tidal damping of eccentricities may 
be a significant stabilizing mechanism over these time 
scales for planets orbiting as close to their star as are 
Kepler-llb and c. Future work should address the limits 
imposed on the eccentricities by the requirement that the 
system remain stable for several Gyr. 

Only two systems became unstable at the nominal 
masses: KOI-284 and KOI-191; KOI-191 is discussed in 
Subsection 15.41 KOI-284 has a pair of candidates with 
periods near 6 days and a period ratio 1.0383 and an 
additional planet at 18 days. Both of the 6-day planets 



would need to have masses about that of the Earth or 
smaller for the system to satisfy Expression ([3]) and thus 
be Hill stable. But the vetting flags are 3 for both of 
these candidates — meaning that they are suspect can- 
didates for other reasons — and the Kepler data display 
significant correlated noise that may be responsible for 
these detections. We expect one or both of these can- 
didates are not planets, or if they are, they are not or- 
biting the same star. This is the only clear example of 
a dynamically-identified likely false positive in the entire 
Kepler multi-candidate population. 

In our investigation, we also were able to correct a 
poorly fit radius of 35 (accompanied by a grazing 
impact parameter in the original model, trying to ac- 
count for a "V" shape) for KOI-1426.03, which resulted 
in too large a nominal mass, driving the system unstable. 
It is remarkable that stability considerations allowed us 
to identify a poorly-conditioned light curve fit. Once the 
correction was made to 13 R^, the nominal mass allowed 
the system to be stable. 

That so few of the systems failed basic stability mea- 
surements implies that there is not a requirement that a 
large number of planets have substantially smaller den- 
sities than given by the simple formula of Equation ([1]) . 
Stability constraints have the potential to place upper 
limits on the masses and densities of candidates in mul- 
tiply transiting systems, assuming that all the candidates 
are planets. 

5. RESONANCES 

The formation and evolution of planetary systems 
can lead to preferred occ upation of resonan t and near- 
resonant configurations (jGoldreichI [T965t iPeald 119761 : 
iMalhotral [1998) . The abundance of resonant planetary 
systems provides constraints on models of planetary for- 
mation and on the magnitude of differential orbital mi- 
gration. 

The distribution of period ratios of multi-transiting 
candidates in the same system is significantly differ- 
ent from the distribution obtained by taking ratios of 
randomly selected pairs of periods from all 408 multi- 
transiting candidates. Most observed planetary pairs are 
neither in nor very near low-order mean motion reso- 
nances; nonetheless, the number of planetary pairs in or 
near mean motion resonances exceeds that of a random 
distribution (Figures [5| - [7| [10| - [T2|). In this section, 
we quantify this clustering and discuss a few particularly 
interesting candidate resonant systems. 

5.1. Resonance Abundance Analysis 

Kepler's multi-planet systems provide a complemen- 
tary data set to exoplanets detected in RV surveys, in 
that RV planets are typically larger and detected on the 
basis of mass rather than radius. Radial velocity mul- 
tiples have the advantage of being minimally biased by 
relative planetary inclinations. Multiply-transiting plan- 
ets have two advantages over RV planetary systems for 
the study of resonances: First, periods are measured to 
very high accuracy, even for systems with moderately 
low SNR, so period ratios can be confidently determined 
without the harmonic ambiguities that affect RV detec- 
tions (jAnglada-Escude et al.M2010l : iDawson fc Fabryckvl 
l2010f ): however, transits can be missed or noise can 
be mistaken for additional transits when the SNR is 
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very low. Second, resonances siRnific antly enhance the 
signal of transit timing variations (Agol e t "all 120051 : 
iHolman fc Murravl[2005t IVeras fc Fbrd 2010), which can 
also be measured accurately, either to demonstrate or 
confidentl y exclude resonance occupation, given enough 
time fe.g.. IHolman et al.llMol ). 

Planets can be librating in a resonance even when they 
have apparent periods t hat are not p e rfectly commensu- 
rate fPcalc 1976; Mar cv et al l 120011: iRivera fc Lissaueii 
120011: iRagozzinc fc Holmanll2010D . Nevertheless, the pe- 
riod ratios of resonant planets will generally be within 
a few percent of commensurability. In the cumulative 
distribution of period ratios (Figure [5]) , we see that few 
planets appear to be directly in the resonance (with the 
important exception of KOI-730), but placing limits on 
resonance occupation requires additional analysis. 

In the Solar System, there is a known excess of satellite 
pairs near mean-motion resonance (MMR) ([Goldreichl 
11965*). and theoretical models of planet formation and 
migration su ggest that this may be the case for exo- 
plane ts fe.g., iMarcv et al.ll2001l : iTerguem fc Papaloizoul 
|2007[) . Can such an excess be seen in the sample of mul- 
tiply transiting Kepler planetary systems? To answer 
this question, we divide the period ratios of the multiply 
transiting systems into "neighborhoods" that surround 
the first and second-order — 1 and j:j — 2) MMRs. 
The boundaries between these neighborhoods are cho- 
sen at the intermediate, third-order MMRs. Thus, all 
candidate systems with period ratios less than 4:1 are in 
some neighborhood. For example, the neighborhood of 
the 3:1 MMR runs between period ratios of 5:2 and 4:1, 
for the 2:1 it runs from 7:4 to 5:2. With this algorithm, 
the first-order MMRs have neighborhoods (represented 
by solid lines in Figure 1101) that are essentially twice as 
large as those for second-order MMRs (dashed lines in 
Figure llOp . The purpose of dividing the period ratio 
distribution into neighborhoods around the major reso- 
nances is to account for the fact that even random period 
ratio s can be close to some ratio of integers (iGoldreichl 
[T965h . 

We define a variable, that is a measure of the dif- 
ference between an observed period ratio and the MMR 
in its neighborhood. In order to treat all neighborhoods 
equally, is scaled to run from —1 to 1 in each neighbor- 
hood. For first-order MMRs, ( is given by 



TABLE 6 

CotJNTS OF PLANET PAIRS FOUND IN EACH RESONANCE NEICHBORHOOD. 



Ci = 3 



1 



V-1 



— Round 



1 



V-1 



(7) 



where V = Po/Pi is the observed period ratio (always 
greater than unity) and "Round" is the standard round- 
ing function that returns the integer nearest to its ar- 
gument. Similarly, for second-order MMRs, is given 

by 

Ca = 3 (:^^ - Round f ) . (8) 



V-l 



V-1 



Table |6] gives the number of planet pairs taken from 
Tables [THl] that are found in each neighborhood. The 
value of C as a function of the period ratio is shown to- 
gether with the data in Figure [TUl We also list neighbor- 
hood abundances considering only adjacent planet pairs. 
This last case primarily removes pairs in the neighbor- 
hood of the 3:1 MMR, but also a nontrivial fraction of 
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Fig. 10. — The value of as a function of period ratio for the 
first order (solid, blue) and second order (dashed, red) MMRs. Also 
shown are the observed period ratios in the Kepler data. 



those neighboring the 2:1 MMR. The fact that signifi- 
cant numbers of non-adjacent pairs are in these neigh- 
borhoods demonstrates how closely packed many of the 
Kepler multi-transiting systems are. 

We calculate C for each planetary period ratio and 
stack the results for all of the resonances in each MMR 
order. The resulting distribution in ^ is compared to a 
randomly-generated sample drawn from a uniform dis- 
tribution in logT', a uniform distribution in V, and from 
a sample constructed by taking the ratios of a random 
set of the observed periods in the multiple candidate sys- 
tems. This third sample has a distribution that is very 
similar to the logarithmic distribution. Figure [TT] shows 
a histogram of the number of systems as a function of C 
for both the first-order and second-order MMRs. We see 
that the results for the second-order MMRs have some 
qualitative similarities to the first-order MMR, though 
the peaks and valleys are less prominent. In the tests 
described below, we also consider the case where all first 
and second-order MMRs are stacked together, as well as 
this stacked combination with the likely resonant sys- 
tems KOIs 730, 191 and 500 (Subsections [El [El and 
15. 5p removed in order to determine the robustness of our 
results to the influence of a few special systems. Finally, 
we study the distribution of adjacent planets only. 

To test whether the observed distributions in are con- 
sistent with those from our various test samples, we took 
the absolute value of ( for each collection of neighbor- 
hoods and used the Kolmogorov-Smirnov test (KS test) 
to determine the probability that the observed sample 
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Fig. 11. — Probability density of the systems as a function of 
f for first order (top) and second order (bottom) mean motion 
resonances. The most common values for ( for first-order reso- 
nances are small and negative, i.e., lie just outside the correspond- 
ing MMR, while one of the least common values, small and positive, 
lie just inside the MMR. No strong trends are observed for second- 
order resonances. The probability density for the logarithmically- 
distributed test sample is a monotonically decreasing function and 
is shown for reference in both plots by white boxes for those values 
of f for which its value exceeds that of the data. 




1^1 

Fig. 12. — Cumulative distribution function (CDF) of the abso- 
lute value of for the combined first-order and second-order res- 
onances (blue). Also shown is the CDF for the logarithmically 
distributed test sample (red). The results of the KS test for these 
two distributions are given in Table [7] along with the KS test re- 
sults for other combinations of resonances and test distributions. 



is drawn from the same distribution as our test sample. 
The results of these tests are shown in Table [T] Figure [T2] 
shows that the probability density function (PDF) for the 
combined distribution of first and second-order MMRs 
differs from the PDF for the logarithmically-distributed 
sample. 

A few notable results of these tests include: (1) the 



distributions observed within the combined first-order 
and second-order neighborhoods are distinctly inconsis- 
tent with any of the trial period distributions — including 
a restricted analysis where the likely resonant systems 
KOIs 191, 500, and 730 are excluded, (2) the distri- 
butions within the neighborhoods surrounding the first- 
order MMRs are also very unlikely to originate from the 
test distributions, (3) there is significant evidence that 
the distributions within the neighborhood surrounding 
the 2:1 MMR alone is inconsistent with the test dis- 
tributions, (4) there is a hint that the distributions in 
the neighborhoods of second-order resonances are not 
consistent with the test distributions, particularly the 
test distribution that is uniform in V. Additional data 
are necessary to produce higher significance results for 
the second-order resonances, and (5) when only adjacent 
planet pairs are considered, all of the test distributions 
are rejected with higher significance. 

Given the distinct differences between our quasi- 
random test distributions and the observed distribution, 
a careful look at the histograms shown in Figure [TT] re- 
veals that the most common location for a pair of plan- 
ets to reside is slightly exterior to the MMR (the planets 
are farther apart than the resonance location), regard- 
less of whether the resonance is first-order or second- 
order. Also, a slight majority (just over 60%) of the 
period ratios have C, values between -1/2 and +1/2, 
while all of the test distributions have roughly 50% be- 
tween these two values in ^ (note that ^ = ±1/2 cor- 
responds to sixth-order resonances in first-order neigh- 
borhoods and twelfth-order resonances in second-order 
neighborhoods — none of which are likely to be strong 
for these planet pairs). There are very few examples 
of systems with planet pairs that lie slightly closer to 
each other than the first-order resonances, and for the 
second-order resonances only the 3:1 MMR has planet 
pairs just interior to it. Finally, there is a hint that the 
planets near the 2:1 MMR have a wider range of orbital 
periods, while those near the 3:2 appear to have shorter 
periods and smaller radii on average. While additional 
data are necessary to claim these statements with high 
confidence, the observations are nonetheless interesting 
and merit serious investigation to determine the mecha- 
nisms that might produce such distributions. 

A preference for resonant period ratios cannot be ex- 
hibited by systems without direct dynamical interaction. 
Multiple star systems with similar orbital periods (i.e., 
non-hierarchical) are not dynamically stable, so dynam- 
ical interaction between objects with period ratios < 5 
indicates low masses, typically in the planetary or brown 
dwarf regime. Thus, the preference for periods near res- 
onances in our dataset is statistically significant evidence 
that most if not all of the candidate (near-)resonant 
systems are actually systems of two or more planets. 
Note, however, that assessing the probability of dynami- 
cal proximity to resonance for the purpose of eliminating 
the false positive hypothesis in any particular system re- 
quires a specific investigation. 

An excess of pairs of planets with separations 
slight ly wider than nominal resonances was predicted 
by Te rquem fc Papaloizou (2007), who pointed out that 
tidal interactions with the star will, under some cir- 
cumstances, break resonances. However, not all of the 
systems near resonance have short-period planets that 
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TABLE 7 

KS TEST P- VALUES FOR DIFFERENT RESONANCE SETS AND FOR TEST DISTRIBUTIONS. 



Resonance set 


logCP) 


V 


All Kepler multis 


2:1 only 


0.00099 


0.00011 


0.00059 


r-j - 1 


0.0012 


0.00021 


0.00044 


j-3 - 2 


0.046 


0.0094 


0.040 


all - 1 and j:j - 2) 


8.2 X 10"^ 


2.7 X 10"^ 


3.9 X lO"'' 


all except KOIs 191, 500, 730 


0.0010 


0.00014 


0.00072 


all adjacent pairs 


9.7 X 10"" 


7.5 X 10"^ 


7.5 X lO"** 



would be strongly affected by tides. Interestingly, KOI- 
730 seems to fiave maintained stability and resonance 
occupation despite the short periods of the planets that 
imply susceptibility of the system to differential tidal evo- 
lution. 

5.2. Frequency of Resonant Systems 

Radial velocity planet searches suggest that roughly 
one-third of the multiple planet systems that have been 
well-characterized by RV observations are near a low- 
order peri od commensurabilit y, with one-sixth near the 
2:1 MMR (jWright et alJl20lH ). Kepler observations find 
that at least ~ 16% of multiple transiting planet candi- 
date systems contain at least one pair of transiting plan- 
ets close to a 2:1 period commensurability (1.83 < V < 
2.18). While it is tempting to consider these results as 
showing a similarity between the Jupiter-mass planets 
that dominate the RV sample and those of the Neptune- 
size planets that form the bulk of the Kepler candidates, 
there are many differences in the criteria used to identify 
these planets and to calculate the period ratios. 

Radial velocity planets are spread over a much wider 
range of period ratio than are the Kepler candidates, and 
the concentration seen in the RV planets is in the region 
very near the 2:1 MMR, whereas the one-sixth number 
for the Kepler candid ates is in the wider " neighborhood" 
that we have defined. iWright et al.l ()201lD 's estimate for 
RV systems near the 2:1 MMR may be lower than the 
actual value for the ensemble of systems that they consid- 
ered because of detection biases. Radial velocity observa- 
tions for systems near the 2:1 MMR have an approximate 
degeneracy with a single planet on an eccentric orbit. 
This degeneracy makes it difficult for RV observations 
to detect a low-mass planet in the interior location of a 
2:1 M MR (|Giuppone et alj 120091: lAnglada-Escude et alj 
[2010h . 

Among RV-discovered systems, only one-quarter of the 
pairs near the 2:1 MMR have an inner planet signifi- 
cantly less massive than the outer planet. The two ex- 
ceptions (GJ 876 c&b, Ara d&b) each benefited from 
a unusually large (> 100) numbers of RV observations. 
Because of this degeneracy, it is possible that many 2:1 
resonant systems with low-mass, interior members have 
been misidentified as eccentric single planets. Thus, the 
true rate of planetary systems near the 2:1 MMR may 
be significantly greater than one-sixth. 

Also, the true fraction of Kepler systems near a 2:1 pe- 
riod commensurability could be significantly greater than 
~ 16%, since not all planets will transit. If we assume 
that planets near the 2:1 MMR are in the low inclina- 
tion regime, then the outer planet should transit about 
axjai ~ 63% of the time, implying that the true rate 
of detectable planets (in size-period space, not account- 



ing for the geometrical limitations of transit photometry) 
near the 2:1 MMR is > 25%. If a significant fraction of 
these systems are not in the low inclination regime, then 
the true rate of pairs of planets near the 2:1 MMR would 
be even larger. 

Based on the Bll catalog, neighboring transiting 
planet candidates tend to have similar radii, and in the 
majority of cases, the outer planet is slightly larger. 
However, smaller planets are more easily detected in 
shorter-period orbits, and a debiased distribution shows 
no preference for the outer planet to be larger (Section 
[2]). The distribution of planetary radii ratios {Rp,o/ Rp.i) 
for neighboring pairs of transiting planet candidates near 
a 2:1 period commensurability (mean motion resonance, 
MMR) is concentrated between 0.8 and 1.25, (see Fig- 
ure [51 dotted curve) . For reasonable assumptions of a 
mass-radius relationship (Equation [T|), this implies that 
most neighboring transiting planet candidate systems 
have masses within ~40% of each other. Fortunately, 
radial velocity surveys can distinguish between pairs of 
planets with similar masses and a single eccentric planet. 
However, those pairs in the tail (~ 20%) for which the 
radii ratio exceeds 1.4 could be difficult to detect with ra- 
dial velocity observations, since for nominal mass-radius 
values, the RV signatures of these pairs would differ from 
that of a single eccentric planet by less than 30% of the 
velocity amplitude of the outer planet (on the decadal 
timescales that are typical of RV observations). Based 
on the planet radii ratio distribution from Kepler, we es- 
timate that the abundance of neighboring planets near 
the 2:1 MMR could be 20% greater than suggested 
by radial velocity surveys due to the difficulty in distin- 
guishing the RV signature of such systems from a single 
eccentric planet. A more precise determination of the 
frequency of resonant or near-resonant systems based on 
Kepler and/or RV surveys is left for future work. 

5.3. KOI-730: A Multiply- Resonant Candidate System 

While few nearly exact mean motion resonances are 
evident in the sample of Kepler planetary candidates, 
one system stands out as exceptional: The periods of 
the four candidates in KOI-730 satisfy the ratio 8:6:4:3 
to ^ 1 part in 1000 or better. This resonant chain is 
potentially the missing link that explains how planets 
that are subject to migration in a gas or planetesimal 
disk can avoid close encounters with each other, being 
brought to a very closely-packed, yet stable, configura- 
tion. Mechanisms that gently ease planetary pairs out of 
such resonant configurations probably account for the ob- 
served preference for pairs to be just wide of resonances, 
as found above. 

This system is difficult to study because of the faintness 
of the target star {Kepler magnitude Kp — 15.34). This 
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faintness, combined with the smah transit depths means 
individual transits are only marginally detected. Indeed, 
this star was not even observed during Q4 because it was 
considered a marginal target until the first of its planet 
candidates was identified. Also, its location within the 
sky is such that during autumn quarters (Q3, Q7, Qll,...) 
its light falls upon Kepler's CCD Module 3, which failed 
early in Q4, so it can only be observed by Kepler about 
70% of the timO The orbital period of planet candidate 
730.03 was initially (Bll) thought be half its presently- 
reported value, putting it in a 1:1 resonance with 730.02. 
Our current best solution discounts every other transit 
of that solution, placing these two planets within the 2:1 
resonance. 

K OI-730.03 shows some sign of weak TTVs through 
Q2 ()Ford et al.l 120111 ) . The remarkably commensurate 
period ratios of these four candidates give us strong con- 
fidence that they all will eventually be confirmed as plan- 
ets orbiting the same star. We are working on a thorough 
analysis of this system, including transit-timing fits and 
long-term stability. 

5.4. KOI-191: A System Protected by Resonance? 

KOI-191 is one of the first five Kepler targets an- 
nounced as showing multiple transiting exoplan et can- 
didates f Borucki et "alll2011a ': Steffe n et al]|2010D . Since 
that time two additional planets have been seen. Two 
planets in the KOI-191 system have a period ratio of 
1.258, and the estimated Hill separation between these 
two planets is A = 2.9 (Table [3]). One of these planets 
has a radius Rp ~ 1 Rjup. Its mass would need to be 
< 17 Mq) to satisfy Expression ([3]). Alternatively, the two 
planets could be locked in 5:4 mean-motion resonance 
that prevents close approaches even if orbits cross. Nep- 
tune and Plu to are locked in 3:2 mean- motion resonance 
of this type (jCohen fc HubbardI I1965D . Relatively fast 
precession, on the order of 100 orbital periods, would be 
needed to maintain the resonance libration, as the ob- 
served period ratio differs from the nominal resonance 
location by 0.6%. We note that the precession time scale 
of the resonant pair of giant planets orbi ting GJ 876 is 
tens of orbital periods (jMarcv et al.ll200"l[) . In short, the 
dynamical stability of this system merits further investi- 
gation. 

5.5. KOI-500: (Near-)Resonant Five-Candidate System 

KOI-500 is a 5-candidate system with periods 0.986779 
d, 3.072166 d, 4.645353 d, 7.053478 d, and 9.521696 d. 
Neighboring pairs of the outer four of these planets all 
have period ratios about 1% greater than those of first- 
order two-body MMRs. Thus they are nearly commen- 
surate, but unless there is an unexpectedly large amount 
(given the small planetary sizes, see Table|4]) of apse pre- 
cession, they are not locked in librating two-body reso- 
nances. Perhaps this is a consequence of diverging tidal 

The Kepler spacecraft rotates four times per orbit to keep the 
sunshade and solar panels oriented properly. Targets are imaged 
on different parts of the focal plane during different orientations. 
The Kepler orbital period is ~372 days, and the data are grouped 
according to the "quarter" year during which observations were 
made. Regular observations were commenced about 65% of the 
way through an orientation; the standard data taken prior to Ke- 
pler's first "roll" are referred to as Ql, and the 10 days of early 
observations done for bright targets are called QO. Subsequent 
quarters are numbered sequentially: Q2, Q3, ... 



evolution, as predicted by Papaloizou & Terquem (2010). 
The three adjacent period spacings between these four 
planets contribute to the statistics of period ratios just 
wide of resonance (small, negative C)- 

Of particular interest, the combinations of mean mo- 
tions 2n2 — 5n3 -|-3n4 1.6 x 10^^ and 2^3 — 6/14 + 4n5 w 
1.3 X 10~^ are so small (smaller than the uncertainties in 
measured values, i.e., consistent with 0) that we suspect 
two 3-body Laplace- like resonances and/or a 4-planet 
resonance may be active and controlling the dynamics 
of this system. 

5.6. KOI-738 and KOI-787: Planets with Period Ratios 

of 9:7 

The pair of planets in KOI-738 have an observed pe- 
riod ratio V = 1.285871, and the pair in KOI-787 have 
V = 1.284008 (Table [J). The ratio 9/7 = 1.285714..., 
so KOI-738 is less than 1 part in 5000 from exact 9:7 
resonance and KOI-787 is less than 1 part in 500 from 
this resonance. The 9:7 resonance is fairly strong if one 
or both planets have a nontrivial e or if their orbits are 
inclined to one another. 

5.7. KOI-657 and KOI-812: Systems where Resonances 

Suggest Missing Planets 

Many Kepler planets likely have resonant companions 
that have not been detected because they do not transit 
or are too small. Such unseen planets may ultimately be 
detected via TTVs that they induce on observed transit- 
ing planets (jFord et al.ll201ll ). Additionally, if they are 
members of resonant chains similar to that seen in KOI- 
730, the period ratio of two observed transiting planets 
in the chain may indicate the presence of missing links. 
If enough chains are ultimately identified, the fraction 
with missing links would yield constraints on the degree 
of non-planarity of multiply-resonant planetary systems. 
We investigate evidence for such chains in this subsec- 
tion. 

The strongest case for a missing link can be made if 
only one unseen planet is needed to complete the chain 
and the links provided by that planet would be much 
stronger than those without it. Neighboring planets in 
KOI-730 have period ratios close to 4:3 or 3:2; the reso- 
nances which correspond to the most significant spikes in 
the distribution of period ratios are 2:1 and 3:2; combina- 
tions of these period ratios can produce two-link chains 
with period ratios of the non- neighboring planets of 16/9, 
2, 9/4, 8/3, 3 or 4. As 2:1 is a first-order resonance and 
3:1 is a second-order resonance, these period ratios can- 
not (without additional evidence from TTVs) be con- 
sidered evidence for intervening planets. In contrast, 
16/9, 9/4 and 8/3 are fifth-order or higher, suggesting 
that they are not important resonances themselves and 
increasing the likelihood that there is an intermediate 
planet. At third-order, a period ratio of 4 corresponds 
to a resonance that is likely to be weak, but cannot be 
as easily dismissed. 

Several Kepler multi-planet systems have observed 
planets with period ratios that are quite close to small 
integer ratios, but with the integers differing by more 
than two, suggesting associations with high-order reso- 
nances that would be quite weak unless the eccentricity 
of one or both bodies is large. However, if one or two 
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additional, as yet unseen, planet (s) orbits at an appro- 
priate intermediate position, then the system could be 
connected via much stronger first-order resonances, as 
is KOI-730 (Section 15.31) . We consider here two specific 
systems in some detail, then list several other KOIs that 
have interesting period ratios, and conclude this section 
with a discussion of situations where the planets may be 
near but not in resonance. 

The observed period ratio of the pair of planets or- 
biting KOI-657 is P = 4.00127 (Table [l]). This ratio 
is within one part in 3000 of the 4:1 resonance. There 
are only 12 neighboring planets with 3.5 < V < 4.5, so 
the a priori probability that any of these period ratios 
would be at least this close to 4 is < 3%. But this res- 
onance is third order, and thus it would require a large 
eccentricity to have significant strength. In contrast, two 
2:1 resonances with an unseen middle planet would be 
far stronger; the Hill separation A between the observed 
pair is 42.5, so there is plenty of room for an intermediate 
planet, even if it is significantly more massive than the 
two that are observed. Note that an alternative chain 
consisting of two intermediate planets in 2:1, 3:2 and 4:3 
resonances, or three intermediate planets involving two 
3:2 commensurabilities plus two 4:3 commensurabilities, 
could also provide linkage via first-order resonances 

An analogous system is KOI-812, with the inner 2 
candidates having a period ratio of 6.005814 (Tabled]). 
There are only 6 neighboring planets with 5.5 < V < 6.5, 
so the a priori probability that any of these period ratios 
would be at least this close to 6 is ~ 7%. The nearby 6:1 
resonance is very weak, but a 3:1 and 2:1 with two inter- 
mediate planets would be significantly stronger, and two 
2:1 resonances plus a 3:2 resonance, requiring three in- 
termediate planets, would be stronger still. The nominal 
A between these two observed candidates is 39, so there 
is enough room for two or three intermediate planets of 
mass comparable to those expected for the observed can- 
didates. 

The candidate systems KOI-82, KOI-117, KOI-124, 
KOI-313, KOI-961 and KOI-1203 also have pairs of plan- 
ets near products of small integer ratios. But these sys- 
tems do not present as strong cases for undetected res- 
onant candidates as the two systems mentioned in the 
title of this subsection because they have one or more of 
the following complicating factors: (i) At least one of the 
intermediate resonances would need to be the 5:4 and/or 
of second order (at most one 5:4 pairing exists among 
observed pairs, described above in Subsection 15.41 and 
it involves a giant planet and the period ratio deviates 
significantly from the ratio of the resonance integers), 
(ii) The estimated Hill separation. A, required if the in- 
termediate planet existed would be small, (iii) There 
is a non-resonant planet candidate within the suggested 
resonant planet chain, (iv) The observed period ratio 
deviates from the small integer ratio by a relatively large 
amount. 

As discussed in Subsection 15. 2[ much of the excess of 
observed near-resonant pairings involves period ratios of 
order 1% larger than small integer ratios. A search for 
chains with this magnitude of deviation from small in- 
teger ratios would require such a broad range of allowed 
period ratios that several pairs would likely be identi- 
fied even if no physical process was favoring such ratios. 
However, such identifications may prove useful in search- 



ing for periods of unobserved planets that might cause 
TTVs in observed candidates. 

6. COPLANARITY OF PLANETARY SYSTEMS 

The orbits of the planets in our Solar System lie close 
to the same plane. The mean inclination of planetary or- 
bits to the invariable plane of the Solar System is z ss 2°; 
if Mercury is excluded, then i « 1.5° and the set of in- 
clinations is similar to a Rayleigh distribution of width 
(Ti « 1°. The strong level of co-planarity has been rec- 
ognized as an important con straint on models of plane t 
formation for over 250 years (|Kantlll755l : ILaolacelllTQfih . 
and our current understanding is that growth within a 
dissipative protoplanetary disk generally yields ci rcular 
orbits and low relative inclinations fe.g.. iSafronoviil969t 
lLissauei]|1993[ ). The interplay between interactions with 
the protoplanetary disk, forcing by distant perturbers, 
and planet-planet scattering, leaves a fingerprint in the 
coupled eccentricity-inclination distribution of exoplan- 
etary systems. High eccentricities are typical for RV- 
detected giant exoplanets whose orbital periods exceed 
one week. Some mechanisms that could produce these ec- 
centric orbits would increase inclinations as well; planet- 
plan et scatterin g is an exa mple of such a pro cess (e.g., 
Chat terjee et a l. 2008; J uric fc Tremainell2008l ). The in- 
clination distribution of typical planetary systems pro- 
vides an important test for planet formation theories. 

The inclination distribution of exoplanets is a fun- 
damental aspect of planetary system dynamics. Yet 
neither transit observations, nor any other technique, 
have directly measured the true mutual inclination be- 
tween planets observed in multiple syste ms, except in un- 
usually fortuitous circumstances (e.g., IWolszczanI 120081 : 
iCorreia et al.ll2010() . This continues to be true in multi- 
transiting systems; even though the inclinations to the 
line of sight of all transiting planets must be small, the 
orbits could be rotated around the line of sight and mutu- 
ally inclined to one another. Indirect constraints, how- 
ever, can be obtained for systems with multiple tran- 
siting planets that make such systems the best probe 
of mutual inclinations around main sequence stars. For 
example, the lack of transit duration var iations of the 
pla nets of Kep ler-9 fH olman et al.l l2010[ ) and Kepler- 
11 ([Lissauer et al. 2011) already place interesting lim- 
its on mutual inclinations to be < 10°, while the in- 
clinations to the line-of-sight of the confirmed planet 
Kepler-lOb and its recently-validated com panion Kepler- 
10c require a mutual in clination > 5° (jBatalha et al.l 
1201 H iFressin et alll2011[ ). The lack of transit timing 
variations can also be used to constrain mutual inclina- 
tions (|Bakos et al.|[2009i) . Occasionally it will be possible 
to measure mutual inclinations from exoplanet mutual 
events (Ragozzine et al., in prep.) or measurements of 
the Rossiter-McLaughlin e ffect of planets in the sam e 
system or by other means (Ragozzine fc Holman| [20To[ l. 

Given the large number of multi-candidate systems 
identified in Bll, we can approach the question of the 
inclination distribution of the population statistically. 
The more mutually inclined a given pair of planets is, 
the smaller the probability that multiple planets tran- 
sit (Figure [131 sec also Ragozzine fc Holman 2010), and 
thus the smaller the fraction of multi-planet systems that 
we would expect to see. Therefore, an investigation of 
coplanarity requires addressing the number of planets ex- 
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Fig. 13. — Multiple transit probability for candidate four and 
five planet systems as a function of mutual inclination from Monte 
Carlo simulation. The curves correspond to individual systems as 
identified in the upper right. For each increment in mutual inclina- 
tion, the multiple transit probability was computed by assigning a 
sky-frame inclination to each planet candidate in the system deter- 
mined by the given mutual inclination with respect to an isotropic 
reference plane (i.e., a random observer's line of sight), and a ran- 
dom nodal angle l|Ragozzine fc Holmanll20lOf l. All orbits were as- 
sumed to be circular. Stellar sizes and masses were taken from 
Table 1 of Bll. A transit was defined as the center of the planet 
passing over any part of the star's disk. At low mutual inclina- 
tion, where the planets are nearly coplanar, the probability of all 
the planets transiting is given by the geometric transit probabil- 
ity of the outermost planet. As the mutual inclination increases, 
the probability quickly decays, as detecting all the planets then 
requires a fortuitous alignment of the observer's line of sight and 
the orbital node of each planet. 

pected per planetary system. Of course, the question of 
planetary multiplicity is also inherently interesting to in- 
vestigate. Here we provide a self-consistent estimate of 
the multiplicity and coplanarity of typical systems us- 
ing a statistical approach involving Monte Carlo gener- 
ation of simulated planetary systems. A similar study, 
undertaken contemporaneously with the second phase of 
our w ork, was recently presented by Tremainc & Dong 
(|2011| ). who take an analytical approach that is in some 
ways more general. 

6.1. Simulated Population Model 
6.1.1. Model Background 

The number of observed transiting planets per star is 
affected by several observational biases. The two most 
notable are that the probability of transiting decreases 
with increasing orbital period and that for a given size 
planet transiting a given star, the duty cycle (and thus 
likelihood of detection) is a decreasing function of orbital 
period, because the fraction of time that a planet spends 
in transit is oc With the current observations, we 

have little to no insight on the non-transiting planets in 
these systems, though eventually investigations of tran- 
sit timing and duration variations (or the lack thereof) 
will be able to place limits on such planets (^Steffe n et al.l 
|2010; Ford et al. 2011). If we knew about the presence of 
every planet, transiting or not, then the problem of deter- 
mining the typical multiplicity would be much easier; in 
this sense well-characterized radial velocity surveys have 



the advantage that they are much less sensitive to the 
inclination dispersion. 

A useful way to address the problem of non-detections 
is to create a forward model of transiting planet detec- 
tions based on a minimal number of assumptions and a 
relatively small number of tunable parameters and then 
to compare the output of this model with the properties 
of the observed systems. The model described here uses 
Kepler's observed frequency of planet multiplicity to esti- 
mate typical values for the true multiplicity of planetary 
systems as well as the mutual inclination between plan- 
ets in the same system. After choosing a distribution for 
the number of planets assigned to a random Kepler star 
(characterized by the parameter Np) and their relative 
inclinations (characterized by <Ti), the model computes 
how many of these planets would be detected, requiring a 
geometrical alignment that leads to transit and a signal- 
to-noise ratio (SNR) large enough to be detected. This 
process is done for a sufficiently large number of simu- 
lated planetary systems (we use 10^ per simulation), and 
then the results are compared to the observed frequency 
of multi-transiting systems. Although imperfections in 
the model and both systematic and statistical uncer- 
tainties in the observations only allow an approximate 
answer, Kepler data enable the first serious attempt at 
using observations to measure the typical values of Np 
(number of planets in a system) and cTi (characteristic 
inclination dispersion) for a sample consisting of a large 
number of planetary systems. 

To isolate the distributions of multiplicity and incli- 
nation, other aspects of the simulated systems (such as 
radius and period distributions) must be consistent with 
the true underlying distributions. This requires some in- 
terpretation of the observed distribution and the iden- 
tification of potential biases. Here, we hold the pe- 
riod and radius distributions fixed and pre-select them 
in a way that the output of the simulation will closely 
match the observed distribution. Our techniques for find- 
ing the debiased dis tribu tions differ fro m those used in 
iHoward et all ()2011[ ) and[^udin| (|2011[ ). but is sufficient 
for the results discussed herein. 

A good example of our methodology is the choice of 
the period distribution. The differences between the ob- 
served and true underlying period distributions are de- 
termined primarily by the geometric bias of observing 
planets in transit, which scales as cx p-2/3. As dis- 
cussed in Bll and confirmed by our investigations, af- 
ter correcting for the geometric bias and discarding very 
short periods (P < 3 days) and large planets {Rp > 6 
R®), the remaining Kepler candidates are generally well 
described by a uniform distribution in log P. The largest 
departure of the observed sample of P > 3 days, Rp < 6 
candidates from a uniform distribution in logP is 
an excess of observed planets near ~20 day orbital peri- 
ods (semi-major axes near 0.1 AU), as discussed in Bll. 
We consider a two-component period distribution: one 
component uniformly distributed in logP and another 
represented in a histogram as a Gaussian in P. The pa- 
rameters of the Gaussian and the relative ratio of these 
two components were minimized by taking the resulting 
period distributions and comparing them to the observed 
distribution with a KS-test. The best fit was a Gaussian 
centered at 6.5 days with a width of 11.3 days (reject- 
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ing without replacement periods below 3 days) where the 
fraction of periods drawn from a uniform distribution in 
log P was 27% and that from the Gaussian in P was 73%. 
These values were determined by maximizing the simi- 
larity of the smooth two-component period distribution 
to the observed periods as measured by a KS-test. This 
two-component distribution gave a KS-statistic of 0.0289; 
for our sample sizes, one would expect a KS-statistic this 
large or larger 61% of the time, even if the two samples 
were drawn from the same distribution, i.e., the analyt- 
ical period distribution is a good and smooth approxi- 
mation to the observed period distribution for planets 
between 3 and 125 days and radii under 6 R®, which 
was the sole purpose of its construction. This analytic 
fit to the observed distribution of transiting planets is 
then debiased for the geometrical probability of align- 
ment by assigning each period a weight proportional to 
p2/3^ This debiased analytical distribution worked well 
in our population simulations, but it should be noted 
that this distribution assumes that the probability of 
each planet transiting is independent of the multiplic- 
ity, which is not strictly true in that systems with higher 
multiplicity must accommodate more planets and tend 
to have slightly longer periods. 

Estimating the true radius distribution of Kepler plan- 
ets also requires a debiasing process. A simple histogram 
(dashed in Figure [H]) of planet radii, Rp, from Bll shows 
that the majority of detected planets have radii around 
~2 Earth radii. Below ~1.5-2 R®, the number of plan- 
ets decreases significantly, which is at least partly due to 
small planets escaping detection due to low SNR. To cor- 
rect for this, we restrict our analysis to planets with radii 
between 1.5 and 6 Earth radii and SNR greater than 20 
(shown in the dotted histogram in Figure \T^ . Each of 
these planets is given a weight of the observed SNR di- 
vided by the radius squared; this weight represents the 
SNR debiased based on planet size. Note that the weight 
was not assigned based on planet-star radius ratio, which 
tends to give extra weight to planets around small stars 
relative to larger stars. A planet radius distribution is 
then taken by randomly drawing from the weighted ob- 
served distribution. This distribution is shown by the red 
histogram in Figure 1141 which represents the underlying 
radius distribution and forms the radius distribution used 
by the simulated population. 

There will be transiting planets drawn from this dis- 
tribution that would not have been detectable by Ke- 
pler in Quarters 0-2 . An inspection of the Bll SNR ratios 
suggests that observational incompleteness sets in below 
an SNR of ~20 (which is calculated from all transits of 
Quarters 0-5 binned together). In a trade-off between 
completeness and increasing the number of planets for 
statistical studies, a less conservative SNR cutoff of 16 is 
taken as the limit for a planet to be "detectable." This 
SNR requirement will be applied to both the observed 
Kepler population and the simulated populations. We 
did not consider the Bll vetting flag in our choice of ac- 
ceptable systems as a large fraction of the multiples have 
not yet been investigated in detail. 

In summary, we are modeling the subset of Kepler can- 
didates that meet the following criteria (based on the 
properties as given in Bll): 

• The planet is detected in Quarter 0, 1, or 2; 
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Fig. 14. — Comparison of the debiased radius distribution used 
in the simulated population (solid red line) to the observed sample. 
The dashed histogram shows the radius distribution of all observed 
planets between 3 and 125 days with no limitations. The dotted 
histogram shows the subset of Kepler planets with SNR > 20 and 
periods between 3 and 125 days. This is the starting point for 
the debiasing process, which increases the relative proportion of 
smaller planets that cannot be detected around every star with 
existing candidates. The black solid histogram shows planets that 
would still be detected with SNR > 20 even if they had radii of 1.5 
Earth radii. The planets in the dotted histogram are weighted by 
SNR divided by the radius squared and the debiased distribution 
(solid red histogram) drawn from this weighted distribution. Here 
we show a random subset of the debiased radius distribution with 
the same number of planets as the dotted histogram to show how 
the relative weighting favors smaller planets. 

• the orbital period is between 3 and 125 days; 

• the radius is between 1.5 and 6 Earth radii; and 

• the SNR (as hsted in Bll for QO-5) is at least 16. 

Applying these cuts to Kepler's observed systems from 
Bll results in 479, 71, 20, 1, and 1 systems with 1, 2, 3, 4, 
and 5 acceptable transiting planets, respectively (Table 
[8]) . Throughout this section, we will be referring only to 
this selected sub-population of qualifying planet candi- 
dates. These 572 qualifying systems are likely to include 
a small fraction of false positives (primarily among the 
single-candidate systems) and our simulation does not 
attempt to account for this. Note that the sixth (outer- 
most known) planet of Kepler-11, planet g, had a transit 
that fell into the data gap between Quarters 1 and 2, so 
this system counts as a quintuple for comparison to the 
simulated population. 

6.1.2. Model Description 

The model starts by taking the full stellar population 
from the 160171 Q2 target stars in the exoplanet pro- 
gram (the list of which is available from MAST), with the 
assumed masse s and rad ii from the Kepler Inpu t Cata- 
log (KIC; Lat ham et al.i 2005: B rown et al.ll2011| ). Stars 
with radii larger than 10 Rq and those with no KIC es- 
timate of Teff or logg are removed, as are those which 
did not have a measured TMCDPP (the temporal me- 
dian of the combined differential photometric precision) 
value, resulting in a list of 153599 stars. The distribution 
of stellar parameters for this subset is for our purposes 
statistically indistinguishable from the distribution in- 
cluding all Kepler stars or the population used by Bll. 
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Fig. 15. — SNR distribution of all candidates with periods from 
3 to 125 days and radii between 1.5 and 6 Earth radii, as taken 
from Table 2 of Bll. The departure from the trend of increasing 
number of planets at smaller SNR is due to observational incom- 
pleteness, which appears to be significant for SNR < 20. As a 
trade-off between completeness and increasing the number of plan- 
ets for statistical studies, an SNR cutoff of 16 (vertical line) is taken 
as the completeness limit for the analysis of Section [6] Not shown 
are 48 candidates with a SNR greater than 100. 

Under the constraints described above, we want to cre- 
ate simulated planetary populations to compare to the 
Kepler observations. In particular, for each simulated 
population, we will compare the number of systems with 
j detectable transiting planets to the observed distribu- 
tion from Kepler. The two main parameters describing 
the simulated population is the average number of plan- 
ets per star, Np, and the inclination dispersion width. 

The first parameter of the simulated population is the 
true multiplicity, Np. Our goal was to employ a mini- 
mal number of parameters, so we focus on distributions 
of planetary multiplicities that can be characterized by 
a single value; an alternative approach of fitting an arbi- 
tra ry distribution of p lanetary multiplicities is presented 
by iTremaine fc Pond (pOll). Three distributions were 
considered for the number of planets: a "uniform" dis- 
tribution, a Poisson distribution, and an exponential dis- 
tribution. In the exponential distribution, the probabil- 
ity of having j + 1 planets is a factor of a smaller than 
the probability of having j planets, with a a free pa- 
rameter. This model had a very poor fit to the data 
(anything that matched the triples to doubles ratio pro- 
duced far too many systems with higher multiplicities), 
for any value of a or di, so we do not consider it further. 

We use Np to refer generally to the average number of 
planets per star and Np^jj if these planets are distributed 
uniformly and Np,\ if they are distributed according to 
a Poisson distribution. The "uniform" distribution pa- 
rameter, Np^u, is the mean number of planets per system, 
which is a fixed number of planets per star if Np^u is an 
integer or an appropriately proportioned mix of the two 
surrounding integers if Np^jj is not an integer. For exam- 
ple, if Np^u — 3.25, 75% of systems would be assigned 3 
planets and 25% would be assigned 4 planets. The Pois- 
son distribution parameter Np_\ assigns the number of 
planets per system based on a Poisson distribution with 
mean Np^\. The frequency of stars with planetary sys- 



tems is considered separately (see below), so we require 
the number of planets to be non-zero; for Np^x ^ 2, the 
distribution differs somewhat from a true Poisson distri- 
bution and Np^x can be less than the average number of 
planets for the ensemble. These distributions were cho- 
sen because they can be represented by a single tunable 
parameter and are reasonable approximations to the true 
expected multiplicity distribution for the kinds of planets 
being considered. 

The second parameter of the simulated population, ct^, 
is the dispersion of the inclinations with respect to a ref- 
erence plane. Planetary inclinations with respect to the 
reference plane in any given system were drawn from a 
Rayleigh distribution. The Rayleigh distribution is the 
appropriate choic e for randomly distribut ed inclination 
angles (see, e.g., iFabrvckv fc WinnI [2009f) and requires 
a single value to define the distribution. The Rayleigh 
distribution is characterized by its width, af, the mean 
value of the Rayleigh distribution with this width is 
(Ji •y/7r/2. Note that this is equivalent to the true mutual 
inclinations between planets in these systems following a 
Rayleigh distribution with width ai\f2. 

Two different inclination distributions were consid- 
ered, which we will call unimodal (ct^.c/) and Rayleigh of 
Rayleighs [oiji). In the unimodal distribution, the value 
for the Rayleigh parameter was taken to be the same for 
all systems in the simulated population, and we denote 
this value as di^u- In reality, the typical inclination dis- 
persion for planetary systems will vary based on the num- 
ber and masses of planets, the infiuence of the protoplan- 
etary disk, and other mechanisms. Thus in a second set 
of simulations, the value of Ui for each system was itself 
drawn from a Rayleigh distribution, so that the distribu- 
tion of the individual simulated planets' inclinations rel- 
ative to their reference planes is a Rayleigh of Rayleighs. 
We denote the width of the Rayleigh distribution for the 
simulated population from which the Rayleigh parame- 
ter for each individual system was drawn by fTi,K. For 
small values of CTj^-R,, most systems are nearly coplanar 
(with a tail of a few systems with significant width); 
for larger values of Oiji^ a small number of systems are 
nearly coplanar, most have inclination widths near ai^ji, 
and a some have quite large inclinations. 

In all cases, nodal longitudes are assigned randomly. 
We assume zero eccentricities as, statistically, the ob- 
served multiplicity rates will not depend on eccentricities, 
assuming that the eccentricities and inclinations are not 
pathologically correlated. Note that one possible covari- 
ance is in systems where multiple planets are eccentric 
and the apses tend to be aligned or anti-aligned. Apsidal 
alignment occurs in the presence of dissipation and this 
could lead to an increased correlation between the prob- 
ability of two planets transiting. This is a second-order 
effect, unless the eccentricities are typically significant 
and non-randomly aligned, and we neglected this effect 
in our calculations. 

After a specific non-zero number of planets is assigned 
to each star, the orbital and physical properties of these 
planets are assigned. First, the orbital period is assigned 
from the debiased analytical population described above. 
Each planetary radius is assigned based on the Kepler ob- 
served radii independently of orbital period, as described 
above. As almost all of the observed multiple planet 
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systems are dynamically stable on long timescales (Sec- 
tion 2]), we impose a proximity constraint that rejects 
simulated multiples that are too closely packed. This 
requires assigning a mass to each planet as well. We fol- 
low the mass-radius relation used above (Equation [1]) , 
i.e., Mp — with masses and radii in units of Earth 

masses and radii. The mutual Hill separation between 
each pair of planets, A, can then be calculated as de- 
scribed in Equations ^ and Although the period 
ratios in observed multiple systems deviate from a ran- 
dom distribution due to clustering near resonance (Sec- 
tion [SJ , this is a minor effect compared to the overall 
geometric bias and the stability constraint that we im- 
pose. 

In practice, the systems are built one planet at a time 
until all the planets assigned to this star (based on Np) 
are selected. The first planet is assigned its period, ra- 
dius, and corresponding mass. Starting with the second 
planet, a new period, radius, and mass are assigned and 
A calculated. If any two planets have A < 3.46 or if 
Ai + Ao < 18 for any three consecutive planets, then 
the planet is rejected on stability grounds and the pro- 
cess repeats until the appropriate number of planets have 
been assigned. Very rarely, over the wide value of Np^\ 
that we are considering, large numbers of planets (more 
than 12) are assigned to a particular star that will not 
fit within the period range and stability criteria; in this 
case, as many planets as can fit in 1000 tries are given to 
this star. 

We can gain additional insight by considering systems 
from radial velocity, where even relatively large inclina- 
tion dispersions do not significantly affect detectability. 
Using the assumed mass-radius relation (Equation[T]) and 
accounting for the average random inclination to the line 
of sight, the size range we are considering corresponds to 
minimum masses of ~ 1.8 — 31.5 M^. In the period range 
that we're considering, the largest number of planets in a 
radial velocity system that satisfy thes e constraints is 5 o f 
the 7 planets found around HD 10180 (jLovis et al.ll2011h . 
Furthermore, while adding more planets (up to ~12) is 
possible theoretically in terms of the conservative dy- 
namical stability requirements laid out above, values of 
Np > 6 are more likely to be long-term unstable in actual 
planetary systems. For this reason, we do not continue 
the calculations to distributions of systems with typical 
numbers of planets larger than Np^jj = 7 and Np^\ = 6. 

At this point, the full orbital and physical characteris- 
tics of the planetary systems have been determined, and 
we can begin to assess observability. First, these plane- 
tary systems are rotated by a rotation matrix in a manner 
equivalent to choosing a random point on a sphere for the 
direction of the normal to the reference plane. In the new 
random orientation, the impact parameter of each planet 
is calculated. If the planets are transiting, the SNR of a 
single transit is calculated assuming a box-shaped tran- 
sit with depth of {Rr,/Ri,Y and d uration calculated from 
Equation (15) in IKippina ()2010l ). corresponding to the 
time the center of planet crosses from one limb of the 
star to the other. A random epoch is assigned and the 
number of transits in a 127-day period (corresponding 
to the duration of Q1-Q2) is calculated. To account for 
the duty cycle of 92%, each of these transits can be inde- 
pendently and randomly lost with 8% probability. The 



total signal is then calculated as the signal for a single 
transit times the square root of the number of observed 
transits. This is compared to the estimated noise over 
the course of the duration of the transit, calculated from 
the TMCDPP values for the randomly chosen star from 
Quarter 2 (see Kepler Data Release Notes available from 
MAST) to determine the simulated SNR. As described 
above, we require that our simulated planets have SNR 
of 9.2 in the simulated Quarters 1 and 2 (corresponding 
to a through-Q5 SNR of 16) in order to be detected. We 
find that for our various simulations, approximately 30- 
50% of simulated transiting planets are "missed" due to 
insufficient SNR. 

The major output of the model is the number of stars 
with J detectable transiting planets, for all j. The next 
step is a statistical comparison between a variety of sim- 
ulated populations (with different choices for Np and ai) 
with the observed distribution from Kepler. An assess- 
ment of whether a particular simulation is rejectable is 
done using an exact (non-parametric) test based on the 
multinomial probability of observing a distribution Oj, 
the number of systems with j planets observed by Kepler^ 
and Ej, the expected number of systems with multiplic- 
ity j generated by scaling to the simulated population. 

First, the two distributions are scaled so that the total 
values of O and E are the same, i.e., '}2,0j = ^ Ej = T, 
where T is the total number of Kepler targets with one 
or more qualifying observed planets (T — 572 in the 
case of 1 < j < 6). Since we will consider various com- 
parisons to observed population (e.g., excluding singly- 
transiting systems (Oi)), we will describe the technique 
for the general case where we are comparing only planet 
counts where jmin < j < jmax- One could imagine 
then generating a very large set of Monte Carlo pop- 
ulations by randomly assigning planets to category j 
with a probability pj = Ej/^ Ej . We have verified 
that the random probability of obtaining a distribution 

~ [^jmi,i1^jmin + lT ■ ■ T^jmax] ^hlS MOUtC CarlO IS 

equivalent to the expected multinomial distribution: 

where T = X^j™"^ The multinomial distribution 
is the generalization of the binomial distribution when 
more than two outcomes are possible. 

Due to the large number of possible populations, even 
the most probable distribution [Ej itself) has a low prob- 
ability of being exactly chosen. To test how likely it 
would be to draw the observed Kepler observations from 
the distribution set by a particular simulated population, 
the significance, S", of the null hypothesis that Oj is in- 
distinguishable from Ej , is determined by calculating the 
sum of all the d istributions of more ex treme probabilities 
than observed ()Read fc Cressidfl988f ): 

M,^. (10) 

Xj-.M^ - <Mo ■ 

If the observed and expected distributions are very simi- 
lar, then there is substantial probability in values that are 
more "extreme" than the observed distribution and the 
significance is high {S ~ 1). If the observed and expected 
distributions are totally different, then the probability in 



Kepler's Planetary Systems 



21 



values more extreme is very small, and the significance 
is very low (S ~ 0). The quantity S* is a measure of how 
good a match the two distributions are Note that this 
method gives the most weight to categories with large 
numbers of objects, but no weight to categories where 
the simulated population produces systems. For this 
reason, the simulated population contains a very large 
number of planetary systems, in order to assign at least 
minimal weight to categories that are rare but possible. 

In practice, we calculate the logarithm (employing Stir- 
ling's approximation of the factorial) of the probabilities 
for a large scaled grid of values centered on Ej mov- 
ing outward until 99.5% of the probability distribution 
is sampled; this is sufficient to calculate whether the hy- 
pothesis can be rejected with very high confidence. Note 
that the expected values from the simulated populations 
Ej are only used in defining the probabilities pj , so the 
scaling of Ej to match the number of observed systems 
does not affect the calculation. We will refer to the fit 
between the Kepler and simulated populations as "ad- 
equate" if the S > 0.05 as, in this case, we cannot re- 
ject with more than 95% confidence the null hypothesis 
that the two populations are drawn from the same dis- 
tribution. Even so, simulations with higher S values are 
considered better fits to the observed distribution. 

When the simulated population predicts systems 
with j transiting planets, yet such systems are observed, 
the multinomial statistic technically goes to 0, i.e., the 
simulated population is completely rejected. In reality, 
we are looking for a single population that fits the major- 
ity of the observed systems. In cases where the simulated 
population contains no quadruply or quintuply transiting 
systems {E4 = i?5 = 0), the two observed Kepler systems 
with four or more planets can be considered outliers to 
the general model. To allow for this situation, the multi- 
nomial statistic is only evaluated over a range of j values 
where the number of expected systems is non-zero (ig- 
noring Kepler observations outside this range, and thus 
limiting the test to fitting systems of lower multiplicity) . 
Since triply transiting systems form a non-negligible frac- 
tion of the population, we do find it reasonable to reject 
outright simulated populations that produce no systems 
with three or more planets. 

In the populations where the number of planets per 
system is drawn from a Poisson distribution, the simu- 
lated population will have significant numbers of systems 
with j > 5 transiting planets, where Kepler observations 
saw no systems. In this case, the multinomial statistic 
can still calculate the probability that the simulated pop- 
ulation matches the Kepler observations. However, the 
calculation of the multinomial statistic increases signifi- 
cantly in difficulty with the number of separate j-planet 
bins in the distribution. For simplicity, we have chosen to 
lump all simulated systems with 6 or more planets into a 
single bin, i.e., comparing 06+ = to £'6+ = Eq+Ei+. . . 
in the multinomial statistic. Throughout, reference to 
J = 6 implies this binning, i.e., we will not distinguish 
between Eq and Eq+. Most of the weight is given to 

^® With some additional minor assumptions, this significance 
will be nearly equivalent to that given by the G-test (similar to 
the categorical test), which relies on calculating the G-statistic 
(G = 20j ln{Oj/Ej)) with no weight given to categories with 
no observed or expected systems. 



the categories with more planets, so the arbitrary choice 
of cutting the distribution off at 6 does not affect our 
conclusions. 

The simulated population assumes that all 10^ Monte 
Carlo stars have planetary systems. Bll report that the 
expected number of planets per star based on the Ke- 
pler observations is about 0.34, though this describes a 
wider range of planets than we are considering here; for 
the limited size range we are considering, the Bll value is 
« 0.2. Note that converting this to the fraction of stars 
with planetary systems (one or more planets) requires 
dividing by the typical multiplicity. 

Our investigation could address the fraction of stars 
with planetary systems, /p, by adding a new free pa- 
rameter and fitting to Oj with < j < 6. We choose 
instead to compute fp by first fitting a scaled popula- 
tion of systems that all contain planets to Oj, j > 0, 
and then calculating the scaled number of systems that 
had planets but did not transit, Eq. Since J2^=o^3 
the number of stars with planetary systems needed to 
match the number of planets observed by Kepler, we 
find fp = Ej=o-^j/160171. That is, we can answer the 
question of the frequency of planetary systems after we 
find simulated populations that match the observed fre- 
quencies of systems where at least one transiting planet 
is detected. Note that for all simulations that have a 
uniform number of planets per star that hosts planets, 
Np,u ^ 2, so in these models the number of planetary 
systems is identical to the number of multi-planet sys- 
tems. When we fit only for observed numbers of multi- 
planet systems, the number of singly-transiting systems 
are underpredicted, and the additional single planet sys- 
tems that would be needed to make up this deficit are not 
included in fp. In contrast, when we use a non-zero Pois- 
son distribution for the true number of planets per star 
hosting at least one planet, then some stars are assumed 
to host one planet, and these are included as planet hosts 
in computing fp. 

6.2. Results and Discussion 

The main results of our simulated populations are 
shown in Table HI We find significant differences in the 
best-fit simulated populations when fitting to Oj over 
the entire range 1 < j < 6 compared to fits that do not 
attempt to match the 479 observed singly-transiting sys- 
tems. We distinguish between the significance with and 
without the singly-transiting systems by defining Si,q 
to be the significance when comparing Oj and Ej for 
1 < j < 6 (i.e., jmin = 1 in Equation fTO]) and 6*2,6 to be 
the significance for 2 < j < 6. Both of these significances 
are listed in Table [8l That table lists the best fit from 
each of the eight possible combinations using Np^u or 
Np,\, cr^,u or (Ji,TC, and jmin = 1 or jmin = 2. Two of the 
extreme simulations (for Np^\ and ct^.k) that are very 
poor fits are also given to show how the Kepler obser- 
vations are able to rule out large inclinations with small 
numbers of planets and coplanar systems with large num- 
bers of planets. Our simulations spanned a much larger 
range than represented in the table, as can be seen in Fig- 
ure [TH which shows a contour plot of S2S that includes 
the best fit to multiple planet systems {Np^u — 3.25, 
(Ji,u = 2°, 5*2,6 = 0.94), and Figure [T^ which shows 
a contour plot of 5*1,6 that includes the best fit includ- 
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ing singly-transiting systems {Np^x — 5.5, ai^n = 15°, 
'S'1,6 = 0.46). Both of these figures demonstrate the trend 
that more planets per star that has planets implies that 
a larger inclination dispersion is required to match the 
Kepler observations. 

When looking at the fits to the entire Kepler sample, 
including singly-transiting systems, the observed distri- 
bution is fit by having large numbers of planets (in the 
specified period and radius range) in most systems and 
only observing a small fraction of them in transit due to 
the large inclination dispersion. It is clear from Table |8] 
that the Rayleigh of Rayleighs inclination dispersion is 
a much better match than the unimodal inclination dis- 
persion to the Kepler sample including singly-transiting 
systems, presumably because it allows more multiples to 
have large enough inclinations that observing only a sin- 
gle planet is probable, while still matching the relatively 
high ratio of triples to doubles {O3/O2 — 0.282). A con- 
tour plot showing the significance Si^e as a function of 
Np^\ and ct^.tj is shown in Figure \W\ We do not include 
constraints from radial velocity surveys, but the higher 
multiplicities (Np > 6) would likely lead to many more 
tightly packed and high multiplicity systems than are 
known. 

Given the small number of degrees of freedom (5 mul- 
tiplicities - 1 scaling - 2 parameters) , it is not surprising 
that there are a large range of adequate fits. As ex- 
pected. Figure [19] shows the degeneracy that increasing 
the number of planets per system requires the inclina- 
tion dispersion to increase in order to provide a good 
match to the Kepler data. Any population within a con- 
tour line is an acceptable fit in that it cannot be rejected 
at 95% confidence that the simulated populations with 
these parameters are significantly different than the ob- 
served Kepler systems. Regions within multiple contours 
provide even better fits. 

Not every star has planets in the size and period range 
that we are considering here; the Sun is an example of a 
star lacking such planets. The number of stars required 
to produce the observed population is given by summing 
over the Ej columns in Table [8] and ranges from approx- 
imately 5000-10000. The mean number of planets per 

star, given as < >= Iqqi-ji ^ in Table |8l range from 
0.182 - 0.196 for simulations that best fit 1 < j < 6. 
When including the required additional population of 
single planets to the simulated populations that best fit 
2 < J < 6, the mean number of planets per star is ex- 
pected to be similar. These values are slightly lower than 
~0.20 planets per star estimated by Bll for the subset of 
planets we are considering here. This small discrepancy 
could be partly explained by completeness issues. We 
chose an SNR =16 cutoff that included more planets at 
the cost of being somewhat incomplete, especially in the 
range 16 < SNR < 20. Correcting for this completeness 
would increase the number of simulated planetary sys- 
tems needed to match the observations, bringing up our 
estimates to that calculated in Bll. 

The simulated population models with unimodal in- 
clination dispersions (cj,;/) do much better when they 
do not attempt to match the number of observed sin- 
gles. The difficulty in simultaneously matching singles 
and multis can be seen qualitatively by comparing the ra- 
tio of doubly-transiting systems to singly-transiting sys- 
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Fig. 16. — Comparison of the period-radius distribution of Ke- 
pler planets to the best-fit simulated population described in Sec- 
tion |6] and Table |8] Black plus signs show the radius-period dis- 
tribution of observed Kepler planets that satisfy the SNR, and 
other criteria for inclusion of the Observed population that is to 
be compared to the simulated populations. The plot only shows 
periods from 3 to 125 days and radii between 1.5 and 6 Earth 
radii, which are the limits discussed in Section \E\ Red X 's show 
the simulated planets from our best-fit population (Np^ij = 3.25 
and cTi^ij = 2.0°), with the number of planets scaled to match the 
observed distribution. One and two-dimensional KS tests show 
that the observe d and simulate d distributions are similar. As dis- 
cussed in Howar detaP l(20TlD and |Youdin (2011), the first 126 
days of Kepler data is not complete below 2 Earth radii and the 
apparent drop ofi^ in planets between 1.5 and 2 Earth radii is due 
to incompleteness effects; the actual population probably continues 
to increase in number for smaller and smaller planets. 

tems (O2/O1 = 0.148) to the ratio of triply-transiting 
systems to doubly-transiting systems (O3/O2 = 0.282). 
These two ratios are quite different, yet the probabil- 
ity of an additional planet transiting depends only on 
the semi-major axis ratio and the mutual inc lination 
(|Ragozzine fc HolmanI [2QToI: iTremaine fc Dongj [2011), 
regardless of the number of planets. We have also found 
that the probability that a planet passed the SNR detec- 
tion threshold is also more-or-less independent of mul- 
tiplicity, i.e., the detectability can be estimated on a 
roughly planet-by-planet basis. (In detail, systems with 
high multiplicities tend to have planets with longer peri- 
ods, which typically lowers the duty cycle, and therefore 
the summed SNR, somewhat.) The simulated popula- 
tions have both the geometric and detectability proba- 
bilities decreasing at an approximately constant rate as 
additional planets are added, generally leading to ratios 
Ej+i/Ej that are nearly constant. 

Unlike in the majority of the simulated populations, 
the observed Ojj^i/Oj ratios are not similar and Ke- 
pler may be seeing an "excess" of singly-transiting sys- 
tems. Note that this cannot be attributed to the ap- 
parent fact that "ho t Jupiters" appear to be singletons 
([Latham et al.l l20TTi ). since we are considering objects 
that are smaller than 6 R0. Some of the excess singly- 
transiting systems could be false positives, but we find 
it unlikely that the majority of the shortcoming can be 
ascribed to bad candidates. An excess of singles and/or 
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Name 










Number 


of systems 


with j 


transiting planets [Ej) 


Significance (5) 


< N > 


fp 















1 


2 


3 


4 


5 


6+ 


'S'l.6 


52,6 






All Kepler 










159210 


791 


115 


45 


8 


1 


1 






~0.34 




Selected Kepler (Oj) 










159719 


479 


71 


20 


1 


1 









~0.2 




Simulated 


Np^u - 5.75 


(Ti 


U = 


10" 


(4694.5) 


469.9 


88.4 


12.2 


1.3 


0.1 


0.0 


0.01 


(<0.01) 


0.189 


0.033 


Simulated 


Np^u - 3.25 


(7 


■V — 


2° 


(5258.3) 


(174.2) 


72.3 


18.9 


1.8 








(<0.01) 


0.94 


0.112 


0.034 


Simulated 


Np,x -4.0 


(7 


,u — 


9° 


(7281.9) 


465.7 


88.6 


15.3 


2.1 


0.3 


0.0 


0.11 


(0.11) 


0.196 


0.049 


Simulated 


Np^x - 1.75 


a 


-.^ — 


1° 


(9481.2) 


(229.4) 


68.9 


18.7 


4.4 


1.0 


0.0 


(<0.01) 


0.51 


0.107 


0.061 


Simulated 


Np^u - 5.75 


(Ti 


TZ — 


10° 


(4826.9) 


472.7 


81.9 


14.2 


2.6 


0.5 


0.06 


0.25 


(0.17) 


0.194 


0.034 


Simulated 


Np.u - 3.25 


(7 


,7^ — 


2° 


(5548.2) 


(199.1) 


72.2 


18.8 


2.0 








(<0.01) 


0.86 


0.119 


0.036 


Simulated 


Np^x - 5.5 


<Ti 


■TZ — 


15° 


(5507.3) 


473.7 


77.9 


15.5 


3.4 


1.0 


0.46 


0.56 


(0.43) 


0.209 


0.038 


Simulated 


Np^x - 2.25 


(7 


,7^ — 


3° 


(8892.9) 


(286.2) 


69.3 


18.0 


4.3 


1.1 


0.3 


(<0.01) 


0.52 


0.130 


0.058 


Simulated 


Np.x - 0.5 




TZ = 


10" 


(76302.0) 


(1960.7) 


87.7 


5.0 


0.3 


0.0 


0.0 


(<0.01) 


<0.01 


0.622 


0.489 


Simulated 


Np^x - 4.5 


a 


,71 = 


0° 


(4134.3) 


(84.8) 


43.1 


24.0 


14.1 


7.9 


7.0 


(<0.01) 


<0.01 


0.121 


0.027 



Note. — Summary of best-fit simulated populations from Section[6] The selected planets have radii between 1.5 and 6 Earth radii, periods between 
3 and 125 days, SNR greater than 16, and show at least one transit in Quarter 0, 1, or 2. The second line shows the number of Kepler systems that 
meet these requirements, and it was this distribution of multiplicities that was compared to the simulated populations. The second column shows 
whether the simulated population in question had a true planetary multiplicity drawn from a uniform distribution (Np^u) with specified mean or a 
non-zero Poisson distribution (Np x)- The third column shows whether the Rayleigh parameter specifying the inclination dispersion from which the 
planets in a given system were drawn was the same for all systems {(Ti^u) or whether the value of this parameter was itself Rayleigh distributed {(Ti^iz), 
as well as the value of ai. See Sect ion 16.1.21 for details. The number of systems with 0, 1, 2, 3, 4, 5, and 6 or more detectable transiting planets is 
shown, along with the significance S (see Equation |10| > of the multinomial statistic when compared to the Kepler observations of 1 < ^' < 6 planets 
(Si^q) and 2 < < 6 planets (52, e)- Higher values of S are better fits, and populations with S values less than 0.05 are not adequate since they can 
be rejected at the 95% confidence level. The simulated population is shown scaled to the observed population (i.e., ^^"^"^^ Oj — Ej) for the 

range of j values implied by the S column that is not in parentheses. Further, the numbers of those classes of systems that are not fit attempted to be 
fit by this group of simulations are shown in parentheses. The best fit for each of the four possible combinations of Np and (Ti is shown for both Si,6 

and 82,6 above the line. Two poorly-fitting simulations are presented below the line. The column labeled < A'" >= — ±^±71 gives the mean number 

of planets per star for the observed and simulated populations and fp = -^^pf^ ' fraction of stars with planetary systems fsee I Youdinlf201 ih . See 
the text for discussion. The Np^u and (7i^u values were tested over the range 2-7 and 0-20°, respectively, as shown in Figure flS] The Np^x s^'cid (Ti^ji 
values were tested over the range 0.5-6 and 0-20°, as shown in Figure [l9| The Np^u and (Ti^-n values were tested over the range 2-6 and 0-10° and 
the Np^x and (Ti^u values were tested over the range 0.5-4.5 and 0-10°, respectively. 
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Fig. 17. — Similar to Figure fT6l but including planets with low 
SNR. Black plus signs again show Kepler planets in the same period 
and radius range, including now all detected planets irrespective of 
SNR. The red X 's show planets in the scaled simulated population 
that were rejected as "undetectable" at present, due to insufScient 
SNR, but which would reach a SNR > 16 in an extended 6-year 
if epJer Mission. Thus, the union of the "+" and "x" distributions 
represents an estimate of what a plot of Kepler planets detected at 
SNR > 16 using 6 years of data will look like. Red dots show simu- 
lated planets that were non-transiting in systems with at least one 
transiting planet. As expected for a relatively thin population, the 
number of non-transiting planets increases significantly at longer 
periods. The distributions of these additional planets depends on 
the simulated population, but the results are qualitatively simi- 
lar for other simulated populations that are adequate fits to the 
observed multiplicities (see Section |6ll. 



very high inchnation multi-planet systems, in the form 
of a population distinct from the one providing the over- 
whelming majority of observed multi-transiting systems, 
is strongly suggested on theoretical grounds for long- 
period planets due to a distinction between systems that 
have had large dynamical scatterings in the past versus 
those systems that have remained relatively calm (e.g., 
iLevison et al.1119981 ). and it may be present in the short- 
period population we are modeling here. 

As we are only attempting to match at most 6 mul- 
tiplicity frequencies Oj, we cannot justify adding a full 
second population to our fits since this would require 4 
highly degenerate parameters. Instead, we investigate 
the fits only to the observed population with two or 
more transiting planets (i.e., 2 < j < 5) using the same 
method described above. Most of the constraint on this 
population is from the ratio of 03/02- It is important to 
note that the comparison statistic S is not modified when 
there are fewer degrees of freedom, so that the larger val- 
ues of S'2,6 compared to Si^e may not imply better fits in 
a statistically significant way. 

The best fit for multis alone (S'2,6 = 0.94) is given 
by a population of systems with Np,u — 3.25 (i.e., 75% 
of systems with 3 planets and 25% of systems with 4 
planets) and with a single Rayleigh inclination width 
of CTi.c/ = 2°, implying a mean inclination of i w 2.5°. 
This population produced an excellent match to the ob- 
servations, though it obviously cannot explain systems 
with 5 or more planets (but see Section [7|) . This pop- 
ulation produced only Ei = 174.2 singly-transiting sys- 
tems. The under-production of singly-transiting systems 
is a trend that is shared with all the best-fit distributions 
with 2 < J < 6, suggesting that as many as two-thirds of 
the observed singles could derive from a different distri- 
bution than the simple population model that matches 
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Fig. 18. — Contour plot of the significance (>S2,6i Equation IIOI I. 
i.e., the probability that a simulated population matches the Ke- 
pler observed population of systems with 2 or more transiting plan- 
ets with radii between 1.5 and 6 Earth radii, periods between 3 and 
125 days, and with SNR of 9.2 or greater (see Section [6] and Ta- 
ble (S]!. The horizontal axis shows the average number of planets 
assigned to all stars in the "uniform" model and the vertical axis 
shows the inclination width as drawn from a unimodal Rayleigh dis- 
tribution di^u for the mutual inclinations in degrees. A plus sign 
marks the best-fit model described in the text and Table |8] Con- 
tours showing S2,6 equal to 0.05, 0.2, 0.4, 0.6, and 0.8 are shown. 
Models outside the thicker 0.05 contour would be rejected as un- 
acceptable fits with an 95% confidence level. Note the clear trend 
that, in order to match the Kep/er observations, the simulated pop- 
ulation must have a higher inclination dispersion if the number of 
planets per star is larger. For this multiplicity distribution, popu- 
lations with mean inclinations greater than i = ai ^7r/2 ft: 5° arc 
adequate, but not as good as fits with smaller inclination disper- 
sions. A contour plot considering the significance 5i_6 for the same 
parameters has lower values in general, and they are shifted such 
that at the same value of Np^jj a higher inclination dispersion (by 
about ~2°) is needed. 

the higher muhiphcities. Figure [TSl gives a contour plot 
of the acceptabihty of simulated populations as a func- 
tion of iVp [/ and Ui^jj for 2 < j < 6. 

Figure also demonstrates the trend that more plan- 
ets require a larger inclination dispersion to match the 
Kepler observations. In this case, the observed distribu- 
tion is fit by having large numbers of planets in most 
systems and only observing a small fraction of them in 
transit due to the larger inclination difference. In all 
the fits we are considering here, this tends to create a 
period distribution of the simulated "observed" planets 
in multiple systems that is much more heavily weighted 
towards periods >40 days than is the observed popula- 
tion. Furthermore, the orbital separation as measured 
by mutual Hill radii, in a plot similar to Figure [9l show 
systems with generally tighter distributions than in the 
Kepler population, as this maximizes the probability that 
nearby planets both transit when mutual inclinations are 
high. Despite these shortcomings, we cannot reject these 
simulated populations with large inclination dispersions 
using the statistical tests that we have performed. We 
leave for future work simulated populations that attempt 
to simultaneously fit the multiplicity, period, and mutual 
Hill radii distributions and which include information of 
planetary multiplicity from radial velocity surveys. 

6.3. Conclusions of Coplanarity Study 
Kepler is providing incredibly powerful insight into the 
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Fig. 19. — Same as Figure FlSl but where the number of planets (in 
the specified range) is drawn from a Poisson distribution (without 
allowing zero planets) with width Np^x and the inclination disper- 
sion is given by a Rayleigh of Rayleighs with width o"; tj. In this 
figure, we show contours of Sifi, i.e., these simulated populations 
were fit to all Kepler transiting planet systems. These populations 
are worse (but still adequat e) fi ts to the observed numbers than 
the model shown in Figure 1181 and show the same trend of in- 
creased inclination dispersion required at higher multiplicity. The 
Kepler observations alone cannot put a strong limit on the true 
multiplicity of planetary systems. 

structure of planetary systems. This simulated popula- 
tion model is the first to explore the multiplicity and in- 
clination distribution of these planetary systems, with re- 
sults summarized in Table[8l Kepler has elucidated a new 
population of planetary systems suspected from RV ob- 
servations: a few percent of stars have multiple similar- 
sized, somewhat- coplanar, 1.5-6 i?0 planets with peri- 
ods between 3 and 125 days. While the distribution of 
Kepler candidates implies that few of these systems con- 
tain Jupiter-size planets, the true multiplicities including 
planets of any size (below the 1.5 limit) will actually 
be higher than the numbers described here. 

There is some evidence that the observed distribution 
may require more diversity than a single homogeneous 
population. To match the numbers of higher multiples 
observed by Kepler requires a population with relatively 
low inclination dispersion; mean inclinations less than ~ 
10° are preferred. If all the systems are drawn from the 
same source population, then the mean inclinations could 
be much higher (though this is not necessary to obtain an 
adequate fit). However, the high planetary multiplicities 
required to obtain adequate fits for high inclinations may 
not be consistent with radial velocity surveys (see also, 
iTremaine fc Dondl201ll ). 

As expected, our simulated populations suggest that 
the majority of the doubly-transiting systems and a sub- 
stantial fraction of the singly-transiting systems are prob- 
ably systems with multiple planets in the size-period in- 
tervals considered for this study (see Figure [T71) . Sys- 
tems with large values of A in Tables [T] and [5] are good 
candidates for systems where additional non-transiting 
planets or undetectably small planets may be missing 
in between observed planets. Similarly, systems with 
small values of A could easily have external planets that 
are not transiting (and possibly even non-transiting in- 
ternal planets). This is consistent with the observation 
that ^15% of candidate planets show significant TTVs 
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and that the frequency of candidates that show TTVs is 
indepen dent of the numb er of candidates observed per 
system ()Ford et al.l 1201 Our simulated populations 
are in general agreement with the corre cted populatio n 
models of Bll. lHcward et aLKpoTTh . and lYoudinI (poTl . 
but suggest that planets may be less abundant than esti- 
mated from radial velocity surveys ([Howard et alJ[2010l : 
IMavor et al. l l20TTh . 

Explaining the formation and evolution of this pop- 
ulation of systems in the context of other types of ap- 
parently distinct systems (e.g., hot Jupiter systems and 
the Solar System) should be a major goal of planet for- 
mation theories. Systems with multiple super-Earth-size 
and Neptune-size planets with periods less than 125 days 
at low relative inclinations suggest interactions with a 
protoplanetary disk th at induced migration while damp- 
ing inclinations fe.g.. [Bltsch fc KlevI I20lil ). Note that 
inclinations can be excited during dis k migration when 
plane ts are captured in resonances (jLee fc Thomme^ 
|2009| ). The scattering process that produces the large 
eccentricities observed for giant planets, which would 
also produce large incli n ations (IChatteriee et al.l 120081 : 
lJuric fc Tre^IIamg l200l ILibert fc TsiganisI 120111) . but 
the distribution of transit durations of Kepler planets 
does not show eviden ce for large eccentricities (see also 
iMoorhead eralllMll) . 

There are good prospects for improving our under- 
standing of the planetary populations observed by Ke- 
pler. Besides continued observations that will increase 
the detectability of planets and find new candidates in 
multiple systems, a more detailed population simulator 
could be developed. For example, although the errors in 
the impact parameters of Kepler planets are usually rel- 
atively large, impact param eters do contain inf ormation 
on relative inclinations fe.g.. lFressin et aI]|201lD that was 
not used in the study presented herein. 

In some systems, the detection of or upper limits to 
TTVs or Transit Duration Variations (TDVs) will con- 
strain the unobserve d component of the relative inclina- 
tions of planets (e.g.. lHolman et al]|2010[ ). Potential exo- 
planet mutual events in Kepler multi-planet systems and 
their constrain ts on mutual incl i nation s will be studied in 
more detail by Ragozzi ne et al.l (j2011| ) . Measurements of 
the Rossiter-McLaughlin effect for planets within multi- 
planet systems would constrain orbital inclinations to 
star's e quator as well as the rnutua l inclination between 
planets (jRagozzine fc Holman|[20Tol) . Each measurement 
of the true mutual inclination in individual systems helps 
fill in the picture for the typical inclination distribution of 
planetary systems. Unfortunately, measurements of true 
mutual inclinations will be observationally challenging 
for most Kepler stars. 

7. HOW RARE ARE PLANETARY SYSTEMS SIMILAR TO 
KEPLER- 11? 

The Kepler-11 {— KOI-157) planetary system has six 
transiting planets, whereas only one other Kepler target 
has even as many as five planetary candidates identi- 
fied to date. Moreover, the period ratio of the two inner 
planets of Kepler-11 is only 1.264, which is the 3'"'' low- 
est ratio among the 238 neighboring pairs of transiting 
if ep/er exoplanets (Figure [SJ. The smallest period ratio, 
1.038, is between two weak (vetting flag 3) candidates in 
the KOI-284 system; this system would be unstable for 



any reasonable planetary masses if indeed both of these 
candidates are actually planets and they orbit about the 
same star, so we suspect that these candidates do not 
represent two planets in orbit about the same star (Sec- 
tionH]). The second smallest period ratio, 1.258, is for a 
pair of candidates in the KOI- 191 system that stability 
considerations strongly suggest are locked in a 5:4 mean 
motion resonance that prevents close approaches (Sub- 
section [531) . Thus, Kepler- lib and Kepler- 11c may well 
have the smallest period ratio of any non-resonant pair of 
planets in the entire candidate list. Additionally, the five 
inner planets of Kepler-11 travel on orbits that lie quite 
close to one another, both in terms of period ratio and 
absolute distance. Kepler-11 clearly does not possess a 
"typical" planetary system, but are systems of this type 
quite rare, or simply somewhat scarce? 

Although its solitary nature makes determining a lower 
bound on planetary systems like Kepler-11 an ill-posed 
question (Kepler could have just gotten lucky), we can 
estimate the number of Kepler- 11 systems that would 
not be seen in transit probabilistically. If all the Kepler- 
11 planets were coplanar (a slight non-coplanarity is re- 
quired by the observed impact parameters), then the 
probability of seeing all six planets transit would be the 
same as seeing Kepler-llg transit, which is 1.2%, suggest- 
ing a simple estimate of 160171 x 0.012, or ~ 1 per 2000 
stars. However, the better-fitting population simulations 
described in Section [6l rarely produced with systems of 5 
or more planets (see Table [Sj), raising the possibility that 
Kepler-11 is not just a natural extension of Kepler mul- 
tiple systems. 

The two candidate multi-planet systems most anal- 
ogous to Kepler-11 are KOI-500 and KOI-707. With 
5 planet candidates, KOI-500 is second in abundance 
to Kepler-11, and four of these candidates are close to 
one another in terms of period ratio (Table [J). The 
target star is significantly smaller than the Sun (~0.74 
R©; Bll), and the orbital periods of the candidates are 
shorter, so if confirmed the planets in this system would 
be the most closely spaced in terms of physical distance 
between orbits of any known system of several planets. 
However, four of the planets orbiting KOI-500 appear 
to be locked in three-planet resonances (Subsection [53]), 
whereas no analogous situation is observed in Kepler-11. 
Unfortunately, KOI-500 is quite faint, so it will be more 
difficult to study than is Kepler-11. The four candidates 
in KOI-707 have periods within the range of the periods 
of the 5 inner planets observed in the Kepler-11 system. 
The multi-resonant KOI- 730 (Subsection 15. 3p is almost 
as closely-packed, but the resonances indicate a qualita- 
tively different dynamical configuration. Period ratios of 
neighboring candidates in the four planet KOI-117 sys- 
tem are all less than two, but nonetheless significantly 
larger than in Kepler-11 and KOI-707. None of the five 
other four-candidate systems is nearly as closely spaced. 
Among the 45 three-planet systems, only two, KOI-156 
and KOI-168, have both neighboring pair period ratios 
less than 1.8. 

Considering only super-Earth-size and Neptune-size 
planets with periods < 125 days, we conclude that 
Kepler-11 appears to be an extreme member of a class of 
very flat, closely-packed, planetary systems. This class 
of systems seems to be signiflcantly less common than 
the classes of planetary systems with a single planet and 
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those that yield most of Kepler's muhi-planet detections 
(Section [S]), but nonetheless accounts for roughly 1% of 
Kepler's targets that have planetary candidates. Kepler- 
11 is a relatively bright star {Kp = 13.7), and two of the 
planets have relatively low SNR ('-^SS), so it may be that 
similar systems that include somewhat smaller planets 
around fainter Kepler targets will be revealed with addi- 
tional data. 

8. CONCLUSIONS 

Analysis of the first four and one-half months of Ke- 
pler data reveals 170 targets with more than one transit- 
ing planet candidate (Bll). While the vast majority of 
these candidates have yet to be validated or confirmed as 
true planets, we expect that the fidelity of this subsample 
of Kepler planet candidates is high (Section [3]) , and the 
small fraction of false positives not to affect the robust- 
ness of our statistical results. Incompleteness of the sam- 
ple due to photometric noise and uncertainties in stellar 
parameters provide additional complications that could 
affect our results in a manner that is difficult to quantify. 
Many of our quantitative findings thus do not have error 
bars associated with them. Nonetheless, multi-transiting 
systems from Kepler provide a large and rich dataset that 
can be used to powerfully test theoretical predictions of 
the formation and evolution of planetary systems. 

Our major conclusions are: 

1. The large number of candidate multiple transit- 
ing planet systems observed by Kepler show that 
nearly coplanar multi-planet systems are common 
in short-period orbits around other stars. This re- 
sult holds for planets in the size range of ^ 1.5 - 
6 R® , but not for giant planets (Figure [T]) . Not 
enough data are yet available to assess its viability 
for Earth-size and smaller worlds, nor for planets 
with orbital periods longer than a few months. 

2. Most multiple planet candidates are neither in 
nor very near mean-motion orbital resonances. 
Nonetheless, such resonances and near resonances 
are clearly more numerous than would be the case 
if period ratios were random. First-order reso- 
nances dominate, but second-order resonances also 
are manifest. There appear to be at least three 
classes of resonance-related relationships evident in 
the data: The most abundant are planet pairs that 
have period ratios from one to a few percent larger 
than those of nearby resonances. Some planet pairs 
as well as chains of three or more planets have or- 
bital periods within one part in 1000 of exact first- 
order resonance ratios. A few nearby pairs of plan- 
ets deviate from exact period ratios by of order 1% 
but appear to be protected from close approaches 
by resonantly librating configurations. We note 
several systems that appear to have particularly 
interesting resonance configurations in Section [5l 

3. Almost all candidate systems survived long-term 
dynamical integrations that assume circular, pla- 
nar orbits and a mass-radius relationship (Equa- 
tion [T]) derived from planets within our Solar Sys- 
tem (Section!?]). This bolsters the evidence for the 
fidelity of the sample. 



4. Taken together, the properties of the observed 
candidate multi-transiting systems suggest that 
the majority are free from false positives. Those 
candidates very near mean-motion resonances are 
most likely to be true planetary systems. Multi- 
transiting systems are extremely valuable for the 
study of planet formation, evolution, and dynam- 
ics (Ragozzine & Holman 2010). 

5. Simulated ensemble populations of planetary sys- 
tems generated to match the observed Kepler mul- 
tiple transiting planet systems tend to underpre- 
dict the observed number of singly transiting plan- 
ets. This provides some evidence for a separate 
population or subpopulation of systems that either 
contains only one detectable planet per star or mul- 
tiple planets with high relative orbital inclinations. 
A third, rarer, group of nearly co-planar densely- 
packed multi-planet systems (Section [7]), also ap- 
pears to be present. Note that these populations 
of planetary systems need not be cleanly separated, 
i.e., there may also be significant numbers of inter- 
mediate systems. 

6. Approximately 3-5% of Kepler target stars have 
multiple planets in the 1.5 R® < i?p < 6 R® and 
3 < P < 125 day range (Section i]). 

7. The inclination dispersion of most multiple planet 
systems in the above size-period range appears to 
have a mean of ^10°, suggesting relatively low mu- 
tual inclinations similar to the Solar System. Many 
singly-transiting systems may come from an addi- 
tional population of systems with lower multiplici- 
ties and/or higher inclination dispersions. 

8. Many Kepler targets with one or two observed 
transiting planet candidate (s) must be multi-planet 
systems where additional planets are present that 
are either not transiting and/or too small to be 
detected (Figure [T7|) . Some of known single can- 
didates already show TTVs (|Ford et al.ll20lH) . and 
we expect that many more will eventually be found 
to have such TTVs and that future Kepler obser- 
vations will reveal large numbers of new candidates 
in existing systems. 

The rich population of multi-transiting systems discov- 
ered by Kepler and reported by Bll have immense value 
both as individuals and collectively for improving our un- 
derstanding of the formation and evolution of planetary 
systems. The Kepler spacecraft is scheduled to continue 
to return data on these multi-planet systems for the re- 
mainder of its mission, and the longer temporal baseline 
afforded by these data will allow for the discovery of more 
planets and more accurate measurements of the planets 
and their interactions. An extension of the mission well 
beyond the nominal 3.5 years would provide substantially 
better data for this class of studies. 
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Discovery mission. Funding for this mission is provided 
by NASA's Science Mission Directorate. The authors 
thank the many people who gave so generously of their 
time to make the Kepler mission a success. D. C. F. 
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